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In  order  to  provide  through-life  support  for  structural  components  fabricated  from  advanced 
composite  materials,  it  is  desirable  to  establish  a  damage  tolerance  methodology.  There  are  many 
difficulties  in  meeting  this  objective,  including  the  multiplicity  of  failure  modes,  the  numerous  types 
of  potentially  significant  defects  which  may  arise  during  manitfacture  or  in  service,  and  the 
sensitivity  to  moisture  arui  temperaatre.  This  paper  focuses  on  the  problem  of  impact  damage  and 
reviews  research  undertaken  <tt  ARL  in  recent  years  related  to  the  analysis  and  testing  of  impact 
damage  composite  laminates.  It  is  shown  that  the  residual  strength  of  impact  damaged  laminates 
reduces  to  an  asymptotic  linut  as  the  size  of  the  damage  increases  and  that  uniaxial  S-N  curves  for 
damaged  lanunates  have  a  generic  shape  with  a  pronounced  threshold  ( endurance)  level. 

Recent  developments  in  the  bonded  repcur  of  impact  damage  are  also  discussed  and  a  simple  repair 
methodology  is  presented.  This  methodology  is  then  illustrated  via  a  laboratory  test  program  and  by 
a  repair  to  a  damaged  FIA-I8  Horizontal  Stabilator. 


DSTO# 


©  COMMONWEALTH  OF  AUSTRALIA  1993 


POSTAL  ADDRESS:  Director,  Aeronautical  Research  Laboratory 

506  Lorimer  Street,  Fishermens  Bend 
Victoria  3207,  Australia. 


Page  Nos. 


1.  INTRODUCTION . 1 

2.  CHARACTERISATION  OF  DELAMINATION  RESISTANCE . 1 

2. 1 .  Mode  I  Characterisation . 3 

2.2.  Mode  n  Characterisation . 4 

2.3.  Discussion  of  Test  Results . 4 

3.  FATIGUE  LIFE  PREDICTION . 5 

3.1.  Test  Results . 5 

3.2.  Prediction  of  Damage  Growth . 6 

3.3.  Generalisation  of  the  Growth  Laws . 7 

3.4.  Remarks . 8 

4.  RESIDUAL  STRENGTH  CALCULATION . 8 

4.1.  Understanding  Static  Failure . 10 

4.1.1.  Finite  Element  Analysis . 13 

4.1.2.  Remarks . 13 

5.  APPLICATION  OF  THE  T*  INTEGRAL  TO  PREDICTING  RESIDUAL 

STRENGTH . 13 

5.1.  The  Finite  Element  Model . 13 

5.2.  Results  and  Discussion . 14 

6.  EFFECT  OF  ANISOTROPIC  PLASTICITY . 1 5 

6.1.  Anisotropic  Plasticity . 16 

6.2.  Result  and  Discussion . 18 

7.  IMPLICATIONS  FOR  HERE  COMPOSITE  REPAIRS . 19 

7.1.  Numerical  Analysis  of  Bonded  Repairs . 20 

7.1.1.  Results  and  Discussion . 20 

7.2.  Experimental  Program . 21 

7.2.1.  Specimen  Fabrication . 21 

7.2.2.  Repair  Fabrication . 21 

7.2.3.  Test  Methodology . 21 

7.2.4.  Results  and  Discussion . 22 

8.  DAMAGED  ASSESSMENT  OF  AN  F/A-18  HORIZONTAL  STABILATOR .  22 

8.1.  Background . 22 

8.2.  Structural  Evaluation . 22 

8.2.1.  Test  Description . 23 

8.2.2.  Results  and  Discussion . 23 

8.3.  The  Finite  Element  Model . 24 

8.3. 1 .  Predicted  Reduction  in  Strain  Due  to  an  External  Doubler ...  24 

8.4.  Repair  Evaluation . 25 

8.5.  Remarks . 25 

9.  CONCLUSIONS  AND  RECOMMENDATIONS . 25 

REFERENCES . 26 

TABLES  1-13 

HGURES  1-16 
DISTRIBUTION 
DOCUMENT  CONTROL  DATA 


-1- 


1.  INTRODUCTION 

Graphite-epoxy  composites  have  many  advantages  for  use  as  aircraft  structural  materials, 
including  their  formability,  high  specific  strength  and  stiffness,  resistance  to  cracking  by 
fatigue  loading  and  their  immunity  to  corrosion.  Thus,  besides  producing  lighter,  mwe 
efficient  structures,  the  use  of  composites  should  provide  more  durable  structures 
compared  with  those  manufactured  from  conventional  aluminium  alloys.  Graphite-epoxy 
composites,  while  having  the  above  advantages,  are  prone  to  a  fairly  wide  range  of 
defects  and  damage  which  may  significantly  reduce  residual  strength  [1-6].  Of  the 
various  types  of  defects,  delaminations  (i.e.  single  or  multiple  internal  cracks  whose 
planes  are  parallel  to  the  surface  of  a  component)  [7-9]  arising  in  service  are  probably  the 
most  insidious  because  they  can  cause  reductions  in  compressive  strength  (up  to  6S%  of 
undamaged  strength  [10,11])  and  are  difficult  to  detect.  Delaminations  may  develop 
during  service  due  to  the  presence  of  excessive  inter-laminar  shear  stresses  or  through- 
the-thickness  tensile  stresses  at  holes,  free  edges,  in  the  region  of  section  changes,  and  in 
bonded  joints.  However,  perhaps  the  most  important  source  of  delaminations  is  impact. 
Such  damage  can  occur  from  dropped  tools  or  from  stones  thrown  up  from  the  runway. 
While  impact  can  cause  a  signiflcant  amount  of  delamination,  often  the  only  external 
indication  is  a  very  slight  surface  indentation  [12,13].  This  type  of  damage  is  frequently 
referred  to  as  'barely  visible  impact  damage'  (BVID).  The  problem  of  BVID  is  of 
particular  concern  because  the  damage  is  unlikely  to  be  discovered  unless  the  region  is 
subjected  to  non-destructive  inspection  (NDI),  such  as  ultrasonic  C-scanning.  However, 
unfortunately,  most  routine  NDI  is  likely  to  be  confined  to  potential  hot  spots  such  as 
critical  joints.  Frequent  full-scale  NDI  is  costly  and  time-consuming. 

It  is  highly  desirable  that  procedures  are  available  so  that  the  possible  occurrence  of 
delamination-type  defects  is  allowed  for  in  the  design  and  certification  of  composite 
aircraft  structures  and  in  the  development  of  approaches  for  through-life  support,  to 
provide  a  rationale  for  setting  inspection  intervals,  particularly  for  highly  stressed  regions 
and  also  provide  repair/reject  criteria  [14,15]. 

This  paper  discusses  the  analysis  and  testing  of  impact  damaged  graphite/epoxy 
laminates,  and  their  ability  to  provide  information  to  establish  a  rationale  for  inspection 
intervals. 
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2.  CHARACTERISATION  OF  DELAMINATION  RESISTANCE 

A  detailed  review  of  the  damage  tolerance  of  impact  damaged  composites  was  given  in 
[1].  The  conclusions  of  this  and  other  subsequent  studies  are  summarised  below. 


•  Compressive  residual  strength  (RS)  following  impact  damage  reduced 

to  less  than  40%  (-3000  to  -4000  pe)  of  the  undamaged  strength 

•  BVID  can  have  a  similar  RS  to  visible  damage 

•  BVID  RS  is  lower  than  that  for  a  drilled  hole  of  the  same  size 
,  Artificial  flaws  have  higher  RS  than  BVID 


•  • 
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•  BVID  produces  an  initial  large  drop  in  the  RS  then  a  gradual 
reduction  to  an  asymptotic  limit 

•  Artificial  disboiuls  also  have  an  initial  large  drop  to  a  constant  RS 

•  Bonded  step  lap  joints  behave  in  the  swie  generic  fashion  as 
monolithic  laminates 

Damage  Logaaon 

•  Work  is  required  on  representative  geometries  including  joints 

•  Work  is  required  on  multiple  impact  locations 

Laminate.  Matrix  State 

•  Moisturised  laminates  have  increased  impact  resistance 

•  The  effect  of  moisture  (hot/wet)  is  to  lower  RS  by  approximately  15 
percent 

Laminate  Stacking  Sequence 

•  Artificial  flaws  are  more  damaging  if  the  outer  plies  are  0° 

•  Dispersed  plies  may  improve  RS 

Component  Size  and  Geometry 

•  Large  components  do  not  always  experience  similar  RS  to  laboratory 
specimens 

Stress-Field 

•  Majority  of  the  work  reported  uses  uniaxial  loading 

•  Consideration  must  be  given  to  representative  bi-axial  loads 

Environment 

•  Relative  RS  increases  under  hot/wet  conditions 

•  Freeze/thaw  cycling  may  not  always  increase  damage 

Prior  Service  History 

•  RS  may  initially  increase  and  then  decrease  under  cycling 

•  Peak  compressive  cyclic  stress  are  the  main  influence 

•  More  information  is  required  on  representative  loading  and  truncation 
of  load  spectra 

Damage  Growth 

•  Artificial  flaws,  when  near  surface,  exhibit  marked  growth 

•  Artificial  flaws  deep  in  laminate  may  not  grow 

•  Impact  damage  usually  grows  only  slightly 

It  should  be  noted  that,  as  shown  in  [1],  experimental  tests  and  analytical  studies  have 
shown  that  for  delaminations  in  (a)  composite  laminates  [33,34,35],  (b)  at  step  lap  joints 
[76,77],  or  (c)  at  mechanically  fastened  composite  joints,  a  stage  is  reached  after  which  a 
significant  increase  in  the  size  of  the  damage  does  not  significantly  reduce  the  residual 
strength.  This  significantly  simplified  the  methodology  for  estimating  critical  damage 
size  and  dramatically  reduces  the  size  of  the  required  static  strength  test  program. 
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Oelamination  in  composites  is  generally  a  mixed  mode  fracture  process  with  inter¬ 
laminar  tensile  and  shear  stresses  at  the  delamination  front.  The  inter-laminar  tensile 
stresses  give  rise  to  a  Mode  I  component  (with  the  associated  energy  released  rate  Gj), 
and  the  inter-laminar  shear  stresses  give  rise  to  a  Mode  II  component  (with  energy  release 
rate  Gu).  When  using  the  energy  release  rate  approach  to  characterise  delamination,  it 
must  be  remembered  that  the  growth  will  often  be  both  mixed-mode  and  non-self-similar. 

In  attempting  to  formulate  suitable  failure  criteria  based  on  the  energy  release  rate 
method,  many  researchers  have  designed  a  variety  of  experiments  for  laminated  plates 
with  the  intention  of  isolating  and  characterising  pure  Mode  I,  pure  Mode  II,  and  mixed- 
naode  failure.  A  selection  of  results  are  discussed  in  the  next  section. 

2.1.  Mode  I  Characterisation 

There  has  been  considerable  success  in  characterising  Mode  I  delaminadons.  The 
methods  most  popularly  used  are  the  Double  Cantilever  Beam  (DCB)  test  and  the  free- 
Edge  Delamination  Tensile  (EDT)  test.  The  results  of  these  tests  on  comparable  material 
systems,  by  various  researchers,  are  briefly  summarised  in  Table  la.  It  can  be  seen  that 
the  values  of  Gu  are  fairly  consistent  within  the  type  of  test  method  employed,  varying 
between  80  and  240  J/m^,  depending  on  the  matrix  material  and  the  nature  of  the  fibre- 
matrix  interface.  In  interpreting  the  values  of  Gu  from  DCB  specimens,  attention  should 
be  given  to  the  surface  morphology  of  the  delamination  area,  since  extensive  fibre 
bridging  [16]  or  intra-ply  cracking  [17]  can  influence  the  value  of  the  energy  release  rate 
obtained  so  that  it  is  higher  than  the  value  that  would  have  been  obtained  if  the 
delamination  was  entirely  confined  within  the  matrix-rich  region  and  the  plane  of  initial 
delamination. 

Table  lb  briefly  describes  some  of  the  experimental  and  data  analysis  techniques  used  by 
researchers  in  estimating  the  value  of  Gu-  For  the  DCB  tests,  thw  methods  of  analysis 
are  commonly  used:  they  are  an  analysis  based  on  beam  theory,  [16,  18-21]  the  area 
under  load  displacement  ciuve  method  [22]  and  a  semi-empirical  compliance  method 
[17].  All  three  methods  give  reasonable  results.  The  edge  delamination  tensile  specimen 
has  received  quite  a  lot  of  attention,  notably  from  O'Brien  [23,24].  Laminated  plate 
theory  and  the  simple  rule  of  mixtures  to  compute  the  stiffness  loss  due  to  delaminations 
have  been  used  to  obtain  an  expression  relating  the  critical  energy  release  rate  Gc  to  the 
critical  strain  Ec  measured  at  the  onset  of  delamination.  The  values  of  Gc  so  obtained  are 
often  only  slightly  higher  than  the  value  obtained  through  DCB  specimens  [22].  This  is 
not  surprising  since,  in  EDT  specimens,  the  delaminations  initiate  at  the  free-edge  where 
high  peel  stresses  are  present,  essentially  resulting  in  predominantly  Mode  I  fracture.  A 
disadvantage  of  EDT-type  tests  is  that  they  often  produce  intra-ply  cracking  as  well  as 
delamination,  resulting  in  complex  fracture  surface  morphologies  rather  than  the  'clean' 
delamination  surfaces  encountered  in  unidirectional  DCB  tests.  An  attempt  to  modify  the 
EDT  specimen  by  introducing  starter  delaminations  in  the  form  of  inserts  has  resulted  in 
considerable  lowering  of  the  value  of  Gk  [25].  Further,  it  has  also  been  shown  that  Gu 
obtained  through  EDT  specimens  is  sensitive  to  differences  in  specimen  widths  [25]  (see 
Table  la). 

The  90®  Ontre  Notch  (CN)  specimen  has  also  been  used  by  a  number  of  researchers 
[22,26].  This  method  relies  on  the  computation  of  the  stress  intensity  factor  Ku  which  is, 
in  turn,  related  to  Gu  through  an  expression  involving  material  constants.  Although  the 
experimental  and  analytical  procedures  of  the  90°  CN  test  are  very  different  from  those 
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of  the  EDT  w  DCB  tests,  the  values  of  Gk  obtained  through  these  methods  agree 
reasonable  well,  perhaps  indicating  that  at  least  f(»^  Mode  I  behaviour,  the  test 
methodologies  are  acceptable. 

2.2.  Mode  II  Characterisation 

The  situation  with  Mode  II  delamination  is  unfortunately  very  different.  There  is  wide 
variability  of  results  for  Gik  in  the  literature  for  comparable  composite  systems.  Table  2 
shows  that  the  variation  of  Guc  obtained  ranges  from  154  to  1042  J/m^.  The  experiments 
designed  to  measure  Guc  are  often  derived  from  tests  to  determine  shear  strength. 
Although  the  Cracked  Lap  Shear  (CLS)  specimen  has  been  used  to  determine  Gno  it  is 
essentially  a  mixed-mode  specimen.  The  proportion  of  Gn  present  in  this  test  has  been 
determined,  through  finite  element  analyses,  to  be  between  75%  and  80%  of  the  total 
system  energy  release  rate  [18,19].  Since  the  Cracked  Lap  Shear  specimen  is  a  mixed¬ 
mode  specimen,  the  value  of  Guc  is  dependent  on  the  failure  criterion  used  [19]  and  the 
accuracy  of  the  finite  element  analysis  in  determining  the  proportions  of  Gj  and  Gn. 

Some  researchers  have  attempted  to  obtain  Guc  more  directly  through  pure  Mode  II  tests. 
Russell  &  Street  [20,27]  have  used  an  edge  delaminated  laminate  specimen  simply  loaded 
as  a  beam  specimen.  Chatterjee's  [32]  tlvee-point  bend  test  was  similar  except  that  two 
delaminations  were  introduced  instead  of  one  and  were  located  in  the  flexural 
compressive  region  of  the  specimen.  Chatteijee,  however,  reported  values  of  Guc  about 
twice  the  value  of  those  obtained  by  Russell  &  Street.  Jurf  &  Pipes  [29]  modified  the 
Arcan  test  fixtures  and  used  them  to  test  a  single  edge  notch  specimen.  Because  of  the 
nature  of  the  geometry  required,  this  specimen  had  to  be  about  90  plies  thick  and  bonded 
to  an  aluminium  fixture.  The  main  advantage  of  the  Arcan  test  was  that  it  enabled  Mode 
I,  Mode  II  and  Mixed-Mode  testing  on  a  single  specimen  type.  However,  there  was  also 
a  tendency  for  some  specimens  to  fail  at  the  fixture-specimen  adhesive  bond,  thus 
rendering  the  test  invalid.  Donaldson  [26]  investigated  Mode  II  failure  of  a  notched  in¬ 
plane  crack  through  a  modified  three-rail  shear  test  arrangement.  It  may  be  argued  that 
since  both  the  growth  of  in-plane  delaminations  and  along-the-flbre  cracks  in 
unidirectional  laminates  are  matrix-dominated,  the  three-rail  shear  test  should  yield  a  Gjk 
that  is  dependent  only  upon  the  matrix  material. 

2.3.  Discussion  of  Test  Results 

Much  work  remains  to  be  done  to  properly  characterise  delaminations  in  composite 
laminates,  particularly  their  Mode  n  and  Mix^-Mode  behaviour.  Various  failure  criteria 
proposed  have  met  with  only  limited  success  in  explaining  mixed-mode  failure.  Much 
evidence  suggests  that  Gnc  is  much  higher  than  Gfe.  In  mixed-mode  fracture,  the 
contributions  of  Mode  I  and  Mode  11  components  to  the  total  energy  release  rate  are 
complicated  and  not  readily  tqtparent  In  particular,  the  individual  quantities  of  Gj  and 
Gu  cannot  necessarily  be  simply  added  together  to  obtain  the  total  energy  release  rate. 
This  is  because,  in  tluee-dimensional  delamination  propagation,  the  growth  of  the  crack 
front  is  often  non-self-similar,  and  furthermore  the  contribution  of  Mode  n  may  become 
significant.  The  analyses  and  experiments  carried  out  to  date  have  concentrate  mainly 
on  the  first  two  modes  of  failure,  although  for  thicker  non-unidirectional  laminates,  the 
contributions  of  Modes  n  and  m  became  increasingly  more  significant  than  that  of  Mode 
I.  Qearly  a  failure  surface  in  three-dimensional  (a  failure  envelope  in  two-dimensional) 
stress  space  needs  to  be  developed.  This  goal  is  presently  not  yet  attainable  since  there 
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does  not  appear  to  be  a  set  of  consistent  data  for  mixed-tnode  fracture  fnHn  different 
investigations.  In  proposing  failure  criteria,  expressions  often  used  are  of  the  form: 


Gk 


(1) 


where  m  and  n  are  constants.  Although  useful  in  curve-fitting  the  experimental  data,  such 
expressions  lack  feasible  physical  interpretations.  Indeed  the  extension  of  such  a  law  to 
three-dimensional  failure,  besides  needing  the  value  Gnio  has  no  physical  backing.  Any 
acceptable  failure  criteria  based  on  energy  release  rates  must  have  a  sound  physical  basis 
and  for  the  case  of  two-dimensional,  self-similar  growth,  equation  (1)  must  reduce  to  the 
form  Gtai  =  Gi  +  Gji.  This  philosophy  has  the  following  implications: 


(1)  A  consistent  test  methodology  must  be  developed  for  both  Guc  and  Gnic,  and  must 
allow  for  the  effect  of  local  constraint. 

(2)  Less  emphasis  should  be  given  to  the  curve  fitting  approach  for  mixed-mode 
fracture. 

(3)  A  detailed  three-dimensional  analysis  should  be  undertaken  for  each  test  method 
in  order  to  determine  the  relative  contributions  due  to  Gi,  Gn.  Gni. 
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3.  FATIGUE  LIFE  PREDICTION 

3.1.  Test  Results  * 

Cyclic  fatigue  tests  play  a  central  role  in  the  assessment  of  fatigue  life  and  the 

specification  of  inspection  intervals  for  metallic  components.  This  has  led  to  similar  tests 

for  composites  with  impact  damage.  It  has  been  found  that  the  S-JV  curves  generated 

from  experimental  data  are  often  flat  over  a  large  range  of  cycles  [33-37].  Tension-  • 

tension  fatigue  loading  is  considered  to  be  less  critical  than  compression-compression  and 

tension-compression  loading  [38  40].  This  paper  will  concentrate  on  the  latter  two 

loading  cases.  For  compression-compression  and  tension-compression  loading,  the 

maximum  residual  compressive  load  as  a  fraction  of  the  static  failure  load  (5)  decreases 

from  1 .0  to  typically  0.6  for  N  in  the  range  1  to  10*,  depending  on  the  initial  damage  size.  • 

The  greatest  rate  of  degradation  appears  to  be  during  the  early  part  of  cyclic  loading,  for 

N  up  to  about  100  cycles.  As  the  number  of  cycles,  and  hence  the  damage  size  is 

increased,  very  little  fiirther  residual  strength  degra^tion  is  observed.  In  fact,  it  has  been 

observed  that  S  =  0-6  may  be  taken  as  the  'fatigue  threshold’  value,  below  which  the  • 

component  may  be  assumed  to  have  'infinite  life.'  Various  researchers  repOTt  threshold 

values  of  between  0.6  and  0.8  [1 1,3?  35,41-44].  5-V  curves  arc  invariably  drawn  from  a 

best  fit  of  limited  experimental  data  often  with  large  scatter  and  their  shapes  vary  from 

straight  slopes  to  some  form  of  curve. 

The  differences  in  fatigue  lives  resulting  from  changes  in  maximum  compressive  stress  or  * 

strain  from  different  researchers  are  summarised  in  Table  3.  From  the  figures  shown  in 

the  table,  it  can  be  seen  that  an  uncertainty  in  the  magnitude  of  the  applied  maximum 

compressive  stress,  for  example  10%,  can  lead  to  a  hundred-fold  difference  in  fatigue 

life.  Of  particular  interest  is  the  considerable  scatter  in  the  data  recorded  by  Ryder  ef  al. 

[35]  and  Potter  [34]  which  show  an  tnder  of  magnitude  of  difference  in  the  fatigue  lives 
of  specimens  at  a  nominally  constant  stress  level. 
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Fcv  service  aircraft  the  in-flight  loads  on  a  component  will  not  be  accurately  known. 
There  is  often  a  significant  variation  between  aircraft,  with  the  result  that  the  stress  levels 
will  always  be  subjected  to  a  degree  of  uncmainty.  We  have  seen  that  a  small  change  in 
stress  can  result  in  a  very  large  change  in  fatigue  life.  This  implies  that  laboratory  testing 
caimot  be  expected  to  provide  accurate  estimates  of  the  fatigue  lives  of  impact  damage 
laminates  and  therefore  cannot  be  expected  to  provide  accurate  information  on  inspection 
intervals.  If  this  is  true  then  the  rationale  for  conducting  such  fatigue  tests  b^omes 
questionable.  Indeed,  there  are  two  natural  consequences  of  this  phenomenon; 

(1)  During  the  designing  of  a  component,  the  stresses  should  be  kept  below  the 
threshold  (endurance)  value,  thus  eliminating  the  possibility  of  a  fatigue  failure 
due  to  BVID. 

(2)  If  an  aircraft  is  already  in  service,  and  the  design  procedure  had  not  included  step 
(1),  then  the  emphasis  for  research  and  development  should  be  moved  towards: 

(a)  Static  strength  testing  to  determine  critical  damage  size  and  loads, 
allowing  for  representative  service  conditions;  and 

(b)  developing  a  rational  reject/repair  criterion  for  the  component  which 
includes  a  methodology  for  determining  inspection  intervals. 


Indeed,  as  a  result  of  this  work  it  is  clear  that  further  research  is  required  in  order  to 
understand  the  physical  reasons  why  the  pronounced  fatigue  threshold  (endurance  limit) 
exists. 


3.2.  Prediction  of  Damage  Growth 

In  order  to  characterise  the  delaminadon  growth  using  fracture  mechanics  techniques, 
several  recent  investigations  (see  Chatterjee  [47];  Mohlin  et  al.  [48];  O’Brien  [49]; 
Ramakumar  &  Whitcomb  [50];  Wilkins  et  al.  [18])  have  adopted  a  growth  law 
relationship  of  the  type  (drzAW)  =  cG"  or  (drz/dAO  =  c(AG)"  where  G  is  the  mixed-mode 
delamination  fracture  energy,  c  and  n  are  constants  determined  from  experiments.  Table 
4  gives  the  value  of  the  exponent  n  as  predicted  in  several  studies  for  different  modes, 
that  is  Mode  I,  Mode  II,  and  Mixed-Mode.  At  this  point  it  is  important  to  note  that  the 
values  of  the  exponent  n  given  in  Table  4  are  obtained  by  fitting  a  suitable  curve  to  the 
experimental  data.  Hence,  in  order  to  get  a  good  correlation  for  predicting  the  growth  of 
a  delamination,  it  is  very  impxntant  to  obtain  experimental  data  which  are  accurate  and 
realistic.  Table  4  reveals  a  large  variation  of  n  for  Mode  I  loading  and  a  relatively  small 
variation  for  Mode  II  loading.  Table  S  gives  the  values  of  n  for  Mode  II  which  are  taken 
from  the  experimental  work  by  Russell  &  Street  [52]. 

Indeed  when  evaluating  the  delamination  fracture  energy  in  different  modes  (that  is  Gi, 
Gn,  and  Gm)  by  analytical  and  experimental  means  there  will  be  inherent  errors  which  are 
somewhat  beyond  control.  As  shown  in  section  2  the  measured  values  for  G|c  vary  from 
80  to  240  J/m2  depending  upon  the  matrix  material  and  the  nature  of  the  fibre-matrix 
interface,  and  the  variation  of  G]c  in  replicate  measurements  appear  to  be  much  greater 
than  ±5%.  For  Gn,  the  situation  is  worse.  The  value  of  Gnc  depends  very  strongly  on  the 
test  procedure  and  the  data  reduction  scheme  adopted  (see  section  2).  Indeed  Table  2 
shows  for  T300/5208  a  variation  in  Gnc  of  between  154  and  876  J/m^.  With  the 
difficulties  involved  in  designing  experimental  methods  for  measuring  pure  Mode  II 
fracture  energy,  it  is  often  very  difficult  to  measure  Gnc  to  an  accuracy  better  than  ±10%. 


We  can  expect  similar  ern»rs  to  occur  in  measurements  of  Gj  and  G\\  in  fatigue  tests.  At 
this  stage  there  have  been  limited  attempts  to  measure  Guk  or  to  determine  the  error 
bounds  on  these  measurements. 


The  effect  that  these  experimental  uncertainties  have  on  growth  will  now  be  examined. 
From  Tables  4  and  S,  take  as  a  first  approximation  nM  for  Mode  II  growth.  Then  the 
10%  uncertainty  in  Gn  gives  rise  to  40%  uncertainly  in  daJdN.  This  situation  is  similar 
in  Mode  I.  Although  the  error  in  Gi  is  smaller,  say  S%,  the  exponent  n  is  larger, 
typically  8  (see  Table  4).  This  makes  the  use  of  these  growth  laws,  to  predict  fatigue  life, 
very  questionable. 

An  accurate  prediction  of  delamination  growth  is  very  much  dependent  on  the  accuracy 
of  measuring  the  values  of  delamination  fracture  energy  during  experiments. 
Furthermore,  the  scatter  in  the  fatigue  test  data  for  composites  gives  rise  to  an  additional 
uncertainty  when  it  comes  to  applying  this  growth  law  in  a  service  situation. 

With  the  prevailing  uncertainty  in  accurately  measuring  the  value  of  Mode  II 
delaniination  fracture  energy,  one  questions  whether  it  is  worth  putting  any  more  effort 
into  conducting  extensive  fatigue  studies  in  order  to  predict  delamination  growth  under 
cyclic  loading. 

The  majority  of  the  tests  to  date  have  concentrated  on  what  was  thought  to  be  pure  Mode 
I  or  Mode  II  growth.  However,  the  growth  of  impact  damage  involves  all  three 
components.  These  conclusions  are  immediately  apparent: 


(a)  A  means  of  characterising  Mode  III  growth  is  required. 

(b)  If  the  uncertainty  in  measuring  Gn  for  a  structural  component  cannot  be 
overcome,  then  alternative  ways  of  characterising  the  growth  of  delamination 
need  to  be  considered. 

3  J.  Generalisation  of  the  Growth  Laws 

There  arc  two  methods  currently  used  to  generalise  the  growth  laws  developed  for  pure 
Mode  I  or  Mode  II  delamination  growth.  These  approaches  can  be  written  as  follows; 


da 

dN 

da 

dN 
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where  G  is  the  total  energy  release  rate  with  components  Gi  and  Gn  whilst  n,  «/,  A.,, 

and  A^are  experimentally  determined  constants.  Equation  (3)  is  empirical  in  nature  and 
can  be  adjusted  until  it  fits  the  data  points  for  a  series  of  Mode  I  and  Mode  II  tests. 
Unfortunately  it  is  entirely  non-physical.  It  infers  that  the  structure  can  subdivide  the 
energy,  i.e.  G,  available  for  crack  growth  into  various  components  Gi  and  Gn  and  then 


use  these  components  in  different  ways.  There  is,  therefore,  a  significant  risk  in  using 
this  approach  in  structures  where  growth  is  neither  pure  Mode  I  new  pure  Mode  II.  Whilst 
equation  (2)  has  a  more  physical  interpretation,  it  must  be  remembered  that  G  is  really  a 
vector  quantity  and  docs  not  equal  Ci  +  Gn  unless  growth  is  self-similar.  We  have  also 
seen  that  the  exponent  n  varies  with  the  nature  of  the  growth  (see  Table  4).  For  mixed¬ 
mode  growth  this  infers  that  n  will  need  to  be  separately  determined.  Indeed  this  cannot 
be  totally  achieved  experimentally  since  for  a  general  problem  G  cannot  be  measured 
from  the  movement  of  the  load  points.  A  combined  experimental  and  finite  element 
investigation  is  therefore  needed  for  any  particular  problem  and  the  need  to  account  for 
irreversible  effects,  and  perhaps  even  rate  effects,  in  the  matrix  material  needs  to  be 
evaluated. 

3.4.  Remarks 

From  a  brief  review  of  studies  coiiducted  to  predicted  delamination  growth  in  composites 
under  cyclic  loading  it  can  be  seen  that  it  is  very  difficult  to  predict  delamination  growth. 
Furthermore,  due  to  inaccuracies  involved  in  measuring  the  values  of  Mode  II  and  Mode 
III  delamination  fracture  energy,  which  in  turn  govern  the  delamination  propagation 
under  cyclic  loading,  current  growth  laws  may  be  of  little  use  under  service  conditions. 

To  aid  in  the  development  of  a  method  for  predicting  the  delaminaiion  growth  in 
composites,  the  experiments  conducted  for  this  purpose  should  be  aimed  at; 

(a)  Minimising  the  errors  involved  in  measuring  the  value  of  delamination  fracture 
energy. 

(b)  A  reduction  in  the  scatter. 

(c)  Development  of  alternative  approaches. 

To  date  evaluation  of  the  fracture  energies  has  generally  been  performed  assuming  linear 
elasticity.  However,  the  matrix  material  is  known  to  be  highly  inelastic  and  to  exhibit  a 
strong  rate  dependence  [4].  In  Mode  I  there  is  a  tendency  for  the  Hbre/matrix  interface  to 
fail  thus  limiting  the  effect  of  these  non-lineapties  [4],  This  is,  in  general,  not  true  in 
Modes  II  and  III,  i.e.  shear  dominated  failure,  and  may  be  one  reason  behind  the  large 
scatter  seen  in  Mode  II  tests.  (Shear  rail  tests  to  determine  the  shear  stress/shear  strain 
response  of  graphite  epoxies  reveals  a  behaviour  similar  to  that  for  the  bulk  matrix 
material  with  extensive  creep,  stress  relaxation,  rate  dependency  and  a  large  inelastic 
strain  to  failure  even  at  room  temperatures.)  It  is  therefore  highly  desirable  that  current 
techniques  be  extended  to  allow  for  these  effects. 

4.  RESIDUAL  STRENGTH  CALCULATION 

There  are  several  techniques  available  to  predict  the  residual  strength  of  an  impact 
damaged  composite.  However,  this  report  will  focus  on  a  generalised  energj'  release  rate 
approach. 

The  energy  release  rate  approach:  For  Mode  I  self-similar  crack  growth  of  a  through 
crack  in  the  absence  of  body  forces,  the  energy  release  rate  G  [54],  can  be  written  as: 
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r*=  J  (Wn,-t,^)ds 


Fm-** 


(4) 


where  is  the  energy  density.  F,  is  a  vanishing  small  closed  path  around  the  tip,  it,  is 
the  conqxMient  of  a  unit  ncnmal  to  the  path  in  the  jc,-  direction,  ti  is  the  traction  vector  on 
the  path  F,  defined  as  ri  =  n,  o^.  and  u,  is  the  component  of  displacement 
Here  W  is  defined  [SS]  as: 

H'  =  ie«C^e,-p,e,(T-T.)-^,e,(Af-Af.)+C,(r.M) 


where  C^u  is  the  stiffness  tenscv,  p  y  and  4  y  are  related  to  the  coefficients  of  thermal  and 

moisture  expansion  and  Ci  is  a  function  of  temperature  and  moisture.  Let  us  now 
consider  the  integral  J,  which  we  will  define  as: 

J.  =  \Wn,+t,^ds  (6) 

J  ax, 

where  r ,  is  the  external  boundary  of  the  body.  There  is  often  a  tendency  to  drop  the 
:<ubscript  s  and  refer  to  J,  as  J.  Using  Green's  theorem  it  follows  that: 


o=y.-  J 


dV 


(7) 


where  the  integral  of  the  second  term  in  equation  (7)  is  over  the  area  between  the  two 
contours.  Thus  J,  will  not  equal  the  energy  release  rate  G  unless  the  area  integral 
vanishes.  At  constant  moisture  and  temperature  J,  may  be  equal  to  G.  However,  in 
general  the  area  integral  will  be  non-zero  and  J„  which  is  measured  experimentally  from 
the  movement  of  the  load  points,  [56]  will  not  equal  G. 


In  service  aircraft  heating  is  often  localized  and  the  moisture  content  varies.  This  will 
produce  a  spatial  variation  in  the  tensor  Qu  with  the  result  that  the  area  integral  will  tend 
to  be  non-zero.  Consequently  when  designing  laboratory  tests  care  should  be  taken  to 
reproduce  the  near  tip  stress,  strain,  temperature  and  moisture  fields  rather  than 
reproducing  the  'global  behaviour.'  This  is  particularly  true  if  the  aim  of  the  test  is  to 
establish  such  quantities  as  the  critical  damage  size  or  the  maximum  permissible  load. 


For  a  three-dimensional  fracture  problem  the  integral  on  the  right-hand  side  of  equation 
(4)  is  no  longer  equal  to  the  energy  release  rate  and  is  referred  to  as  t*  [54,56]. 


The  application  of  this  iqiproach  will  be  considered  in  sectimi  6.  It  must  be  noted  that,  as 
mentioned  in  section  2,  experimental  tests  and  analytical  studies  have  shown  that  for 
delaminations  in  (a)  composite  laminates,  (b)  at  step  lap  joints,  or  (c)  at  mecharucally 
fastened  composite  joints,  a  stage  is  reached  after  which  a  significant  increase  in  the  size 
of  the  damage  does  not  signincantly  reduce  the  residual  strength.  This  significandy 
simplified  the  methodology  for  estimating  critical  damage  size  and  dramatically  reduces 
the  size  of  the  required  static  strength  test  program. 

Unlike  the  situation  for  fatigue  life  estimation,  fracture  mechanics  based  tools  provide  a 
useful  tool  in  both  understanding  and  predicting  the  static  failure  process.  To  further 
illustrate  this,  let  us  now  consider  the  analytical  solution  for  a  circular  disbond  in  an 
isotrt^c  plate  under  remote  hydrostatic  compression  [75]. 
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4.1.  Understanding  Static  Failure 

When  a  composite  material  is  subjected  to  impact,  considerable  damage  can  be  caused 
since  the  cohesive  strength  of  the  plate  through  the  thickness  is  very  low  and  this  in  turn 
can  lead  to  a  deterioration  in  the  strength  of  the  plate.  This  failure  pnx:ess  can  be  divided 
into  two  phases.  In  the  first  phase,  the  plate  is  impacted  and  the  resulting  response  causes 
interlaminar  separation.  The  size  of  this  damaged  area  is  a  function  of  the  impact 
parameters,  impact  energy,  plate  material,  lay-up,  etc.  In  the  second  phase,  the  damaged 
area  spreads  to  the  undamaged  area  through  a  combination  of  laminate  buckling  and 
further  delamination. 


In  order  fix’  growth  to  occur,  work  of  rupture  is  required.  This  energy  n:ust  be  available 
from  either  the  work  done  by  the  compression  forces  (denoted  by  W)  during  the  growth 
of  the  delaminadon,  and/or  from  the  decrease  in  the  strain  energy  U  of  the  system. 
Further  delaminadon  growth  after  the  occurrence  of  buckling  will  depend  on  the 
magnitude  of  the  fracture  energy,  r^,  which  is  defined  as  the  energy  requir^  to  produce 
a  unit  of  new  delaminadon  area.  If  the  energy  release  rate  G  at  a  given  instance  is  greater 
than  Tg,  then  further  delar-  Ion  growth  will  occur.  To  understand  this  process, 
consider  a  layered,  isotropic  circular  plate  containing  a  thin-film  delaminadon  of  radius  a 
as  shown  in  Figure  1.  The  plate  is  subject  to  an  in-plane  compression  undl  the  load 
carried  by  the  delaminated  region  reaches  a  cridcal  value,  at  which  buckling  of  this  part 
of  the  structure  occurs,  leading  to:  (i)  the  onset  of  transverse  deflecdons  of  the 
delaminated  region;  and  (ii)  reducing  the  amount  of  further  load  the  delaminadon  can 
support.  This  model  is  defined  as  "Specimen  A"  -  see  Figure  2(a). 


The  delaminated  region  will  buckle  when  the  stress  it  carries  o’  reaches  the  critical 
buckling  stress,  o'*  ,  for  a  clamped  plate  of  radius  a  and  thickness  t  [57],  i.e.,  when; 

o'  =  a'*  =  y^D/a^t  (8) 

where  IT  =  3.832  and  D=£/^/l2(l-v^).  At  this  point,  the  externally  applied  stress  is; 

a  =  o*  =  a'*  (9) 

Prior  to  local  buckling,  the  radial  displacement  at  r=b  is: 

u(b)  =  «,(a)  =  (l-v)o6/£.  (10) 

An  analysis  of  the  system  energy  release  rate  can  be  obtained  [58]  by  using 

2na  da 

At  the  initiation  of  local  buckling,  the  applied  stress  o*  and  the  displacement  at  the 
boundary  of  the  plate  is  u*(fr),  where  it  can  be  shown  [59]  that: 

u*(fe)=(l-v)o*6/£  (12) 

After  buckling,  the  stiffness  of  the  delaminadon  decreases  from  £'=£  to  £'=3i  £,  where 
\  is  given  by  [58]  as: 

X=(l-H.207(l-)-v))-'  (13) 


Because  the  delaminadon  loses  some  of  its  stiffness  due  to  buckling,  this  causes  the 
compliance  of  the  disk  re[0,<i]  to  increase.  This  can  be  simply  modelled  as  a  plate  under 
ctxnpression,  with  an  inner  section  having  Young's  modulus  £|  (<£)  and  an  outer  section 
of  modulus  £  as  shown  in  Figure  2,  "Specimen  B.” 
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For  these  two  specimens  to  be  equivalent,  we  require  the  displacements  at  r=a  to  be 
equal.  i.e.: 

UJ^(a)=u,{a),  (H) 

where  UA(a)  is  the  displacement  at  r-a  for  q)ecimen  A  and  from  [59], 

u^(a)  =  [(1  -  v)o,ti  /  E][\  -  (f  /  rXl  -  X)]"'  (15) 

and,  f(CM’ specimen  B 

Ug{a)  =  il-v)a^a/ E^.  (16) 

Here  ua(o)  was  obtained  from  the  requirement  that 

o,  =  (r/r)o.,-K(r-0/r)Oa  (i7) 

and 

Uai(o)  =  u^(ci)  (18) 

where 

UM(o)  =  ii-v)ana/  E'  (19) 

and 

u^2(a)  =  (l-v)Caa/ E.  (20) 

In  the  above,  UAt(a)  and  o,,  are  the  displacement  and  stress  at  r=a  of  region  1  in 
specimen  A  where  regions  1  and  2  are  defined  in  Figure  2(a).  From  equation  (14)  we 
find  that: 

E,=XE  =  (l-(t/T)(\-X))E.  (21) 

Fot  specimen  B,  the  displacement  at  any  point  €  [a,b\  is  given  by  [59]: 
u(r)  =  -ia^b^lb^-a^) 

x[(-l/rXo-a;)(l-t-v)  (22) 

+r{aJb^-a/a^){l-v)]/E 

To  find  the  displacement  at  the  boundary  n=b  of  specimen  B,  we  require  the  inner 
di^iacement  (at  r=a)  between  regions  1  and  2  to  be  continuous, 

M„(a)  =  aj.(fl),  (23) 

where  Uai  and  ub2  are  given  by  equations  (16)  and  (22).  Solving  equation  (23)  3delds: 

o,  =  o[l-fccr\  (24) 

where 

tr  =  (l-X)/X,  c  =  -(l-v)(b*-a*)/b^  (25a,b) 

2 

Substituting  for  Oj  into  equation  (22)  yields: 

M(6)  =  M,(o)[l  +  e],  (26) 


where 


(27) 


e  =  2ea*  (26^  +  e(l  -  v){b^  -  a" ))-' . 

Thus  the  diqilacement  undergone  since  buckling  is: 

m,(6)  =  (1-vXo-o*)6(1  +  ©]/£  (28) 

and  the  total  displacement  is  given  by: 

u{b)  =  u*{b)-^Uf(b) 

=  «.(o)[l+(l-o*/a)e].  (29) 

Because  of  buckling,  not  only  does  the  compliance  of  the  inner  disk  r  e  [O.a]  increase, 
but  also  that  of  the  whole  specimen.  As  a  result  of  this,  the  load-displacement  curve  is  no 
longer  linear,  but  in  fxx  bi-linear.  (A  typical  load-displacement  curve  is  shown  in  Figtue 
3.) 

With  these  assumptions  it  follows  that: 

G  =  (l-X)(l-v)^(o^-o*^-^m"(o  -a*)^)/E(}  +  ec),  (30) 

where 

m=e(l  -  v)a^  (2h'  +e(b^  -a^)(l  -  v))’*  (31) 

Now  ec  and  m  are  of  the  orders  (r/T)  and  {tcP-fTtP-)  respectively,  thus  for  {tfT)  small, 
equation  (30)  can  be  approximated  by: 

G  =  (l-X)(l-v)r(a^-o*')/£,  (32) 

Plots  of  G  as  a  function  of  the  delamination  radius  a  and  the  load  P  for  a  specimen  of 
dimensions  fr=75  mm,  r=3.1  mm,  ^=0.5  mm  and  modulus  £==68.9  GPa  and  v=0.3  are 
shown  in  Figures  4  and  5. 

In  order  to  generalise  the  energy  release  rate  for  arbitrary  parameters  and  fracture  energy 


r„,  we  introduce  the  following  normalizations: 

r,=r,£/(l-X)(l-v)/,  (33) 

o  =  a/r.''',  (34) 

ifD.  (35) 

With  this  notation  equation  (32)  can  be  rewritten  as: 

=  (36) 


This  solution  reveals  that  as  the  size  of  the  delamination  increases,  the  energy  release 
rate,  and  hence  the  residual  strength,  tends  to  an  asymptotic  value.  This  phenomena  plays 
a  central  role  in  the  damage  tolerance  assessment  of  impact  damage,  and  is  frequently 
observed  in  residual  strength  tests  on  impact  damaged  structural  components,  see  sections 
2  and  4. 


4.1.1.  Finite  Element  Analysis 

Befnxv  perfuming  a  detailed  finite  element  analysis  of  an  impact  damaged  structure  it  is 
first  necessary  to  establish  that  the  analysis  technique  is  capable  of  reproducing  known 
analytical  solutions.  The  problem  described  above  represents  one  such  potential 
solution.  Let  us  consider  the  specific  case  when  the  structure  is  a  symmetrically 
delaminated  qiecimen  with  thickness  IT,  having  dimensions  h=7S  mm.  'f=3.1  mm,  r=O.S 
mm,  and  modulus  £=^.9  GPa  and  v=0.3. 

Due  to  symmetry  it  was  only  necessary  to  model  half  of  the  specimen.  The  finite  element 
mesh  is  shown  in  Figure  6.  This  mesh  was  generated  using  four-noded  axisymmetric 
elements.  Due  to  concern  about  the  adequacy  of  four-noded  elements  to  model  the 
bending  that  occurs  when  the  delamination  bucldes,  another  mesh  was  created,  similar  to 
the  one  above,  except  that  the  delaminated  layer  is  now  composed  of  eight-noded 
elements. 

Four  delamination  radii  a=15,  20,  35  and  45  mm  were  considered,  and  analyses  were 
conducted  to  determine: 

(i)  the  energy  release  rate  as  a  function  of  applied  load; 

(ii)  the  energy  release  rate  as  a  function  of  delamination  radius. 

In  each  case  the  stress  intensity  factor  and  hence  G  was  calculated  from  the  resultant 
opening  of  the  faces  of  the  crack. 

The  resultant  values  of  G  are  shown  in  Figures  7  and  8,  as  a  function  of  applied  load  and 
radius  of  the  delaminadon,  together  wdth  the  corresponding  analytical  values.  From  these 
figures  it  can  be  seen  that  the  closed  form  analytictd  solution  for  G  agrees  quite  well  with 
the  finite  element  solution,  obtained  using  eight  noded  elements. 

4.12.  Remarks 

A  simple  closed  form  solution  for  the  energy  release  rate  G  has  been  obtained  for  a 
circular  delaminatioa  in  a  circular  plate  and  was  verified  by  means  of  finite  element 
analysis.  This  solution  reveals  the  dependence  of  G  on  the  geometry  and  location  of  the 
delamination.  It  may  also  be  useful  as  a  "benchmarit"  test  for  evaluating  the  ability  of 
various  finite  element  packages  to  model  delamination  damage. 

It  is  important  to  note  that  as  the  size  of  the  delamination  increases,  G  tends  to  an 
asymptotic  limit  This  is  consistent  with  observed  experimental  results  for  which  the 
residual  strength  asymptotes  to  a  lower  bound  as  the  size  of  the  impact  damage  increases. 


5.  APPLICATION  OF  THE  T*  INTEGRAL  TO  PREDICTING  RESIDUAL 
STRENGTH 

5.1.  The  Finite  Element  Model 

Let  us  now  turn  to  the  problem  of  a  laminated  composite  structure  containing  an  internal 
delamination.  Details  of  the  finite  element  mesh  us^  in  this  work  are  shown  in  Figure  9. 
It  is  a  model  of  an  impact  damaged  laminate  with  a  fastener  hole  under  compression. 
The  dimensions  of  the  model  are  the  same  as  those  used  in  the  experimental  work  of  [35]. 
The  specimen  tested  as  a  [0/45/02/-45/02/45A)2/'45A)]t  T300/5208  graphite-epoxy 
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Uuninate  and  contained  a  centrally  located  hole  of  9.5  mm  diameter,  surrounded  by 
delamination  damage  due  to  impact  and  poor  drilling. 

The  elements  used  ate  mostly  twenty-noded  iso-parametric  elements  with  directionally 
reduced  integration  and  2x2x3  Gaussian  quadrature  points,  with  the  3  points  being  taken 
through  the  ply  thickness.  Detailed  description  of  the  reduced  integration  scheme  can  be 
found  in  [60].  The  crack  tip  elements  along  the  circular  delamination  are  fifteen-noded 
iso-parametric  wedge  elements.  The  effect  of  moving  the  relevant  mid-side  nodes  of 
these  elements  to  the  quarter  point  positions  to  simulate  the  desired  singularity  is 
investigated.  Solutions  with  and  without  quarter  point  shifts  are  compared. 

The  initial  damage  around  the  fastener  hole  from  [35]  is  modelled  as  a  circular 
delamination  of  radius  13.75  mm  between  the  second  and  third  plies  (i.e.  between  the  45° 
and  0°  plies).  This  is  an  approximate  simulation  of  the  initial  damage  as  shown  in  Figure 
10(a).  It  can  be  seen  hrom  the  ultrasonic  C-scan  that  the  initial  delamination  is  nearly 
circular. 

The  two  plies  above  and  below  the  delaminadon  are  modelled  separately  with  ordinary 
three  dimensional  elements  while  the  remaining  20  plies  are  modelled  with  super¬ 
elements  with  displacements  varying  quadratically  in  the  local  iso-parametric  coordinate 
system.  The  material  properties  used  are  those  of  Narmco-Thomel  T300/5208  and  are 
tabulated  in  Table  6.  The  effect  of  changing  the  value  of  the  shear  modulus  C/23  is 
investigated.  In  the  first  case,  the  value  of  C/23  is  taken  as  equal  to  that  of  G12  and  C13.  In 
the  second  case,  G23  is  estimated  via  the  approximation: 

GJ,=£J^/2(l-^Va)  (37) 

which  was  also  used  in  [61].  The  properties  of  the  matrix  material  are  given  in  Table  6. 

It  is  important  that  in  the  finite  element  model,  the  faces  of  the  delamination  are 
prevented  from  overlapping.  Otherwise,  non-physical  solutions  may  be  obtained.  By 
examining  the  solutions  of  the  displacements,  it  is  found  that  some  parts  of  the 
delaminated  faces  have  overlapped.  Thus  a  series  of  constraint  equations  are  applied  to 
tqtpropriate  nodes  to  simulate  local  closure.  The  effect  of  correctly  restraining  the 
delamination  faces  is  also  studied,  with  three  sets  of  solutions  involving:  (i)  incorrect 
restraints,  (ii)  correct  restraints,  and  (iii)  no  applied  restraints.  A  compressive  strain  of 
0.(X)6  is  applied  at  the  ends  of  the  model.  This  value  of  strain  represents  the 
experiment^y  measured,  far-field  failure  strain  for  this  specimen,  see  [35]. 

5.2.  Results  and  Discussion 

An  interesting  characteristic  of  the  distribution  of  T*,  in  this  particular  problem,  is  the 
existence  of  three  peaks  along  the  circular  delamination  front,  see  Figiue  11.  The  T* 
distribution  has  the  largest  peak  located  at  0=60°,  where  0  is  defined  in  Figure  9.  The 
value  of  T*^  for  the  correct  restraint  condition  is  70.8  Jm*^.  This  value  compares 
favourably  with  the  Gic  value  of  76  Jm"^  for  neat  5208  Narmco  resin  obtained  by  Bascom 
et  al.  [62].  Double  cantilever  Mode  I  tests  provide  values  of  78-88  Jm'*  for  the  Gk  of  the 
T300/5208  material  system  [17,18]. 

It  must  be  mentioned  that  the  value  of  T^max  is  very  sensitive  to  the  restraint  condition 
iq>plied  to  the  delamination  faces  in  order  to  prevent  overlapping  of  the  surfaces[48]. 
Incorrect  conditions  result  in  an  1 1%  error  in  T*. 
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The  influence  of  shear  nxxiuli  is  also  shown  to  be  significant,  particularly  in  the 
distribution  of  T*.  This  effect  is  imptniant  since  the  value  of  G23  us^  in  the  numerical 
simulation  affects  the  calculation  of  interlaminar  shear  stresses  crucial  in  the 
determinatitMi  of  delamination  behaviour.  In  Figure  12  it  can  be  seen  that  although  the 
basic  shape  of  the  distribution  of  T*  remains  quite  the  same,  the  relative  and  absolute 
magnitude  of  the  peaks  are  greatly  affected.  The  two  cases  considered  are  perhaps 
extreme  cases  since  the  true  ot  effective  value  of  G23  may  lie  somewhere  in  between  the 
values  used  in  cases  A  and  B. 

When  quarter  point  shifts  of  the  crack  tip  nodes  are  not  performed  the  distribution  of  T* 
is  changed  remarkably  (see  Figure  13).  In  this  case  the  analysis  failed  to  detect  the 
secondary  peak  at  0=45°.  These  results  clearly  illustrate  the  importance  of  correctly 
modelling  the  stress  singularity  and  the  need  to  prevent  overlapping  of  the  delamination 
faces. 

Figure  10b  shows  the  c-scan  of  the  final  shape  and  size  of  the  damage  after  20,000 
cycles.  It  can  be  seen  that  the  points  of  maximum  growth  coincide  with  the  local  maxima 
of  T*.  Furthermore,  as  mentioned  earlier,  the  predicted  value  of  T*„^  70.8  Jm'^,  at  the 
failure  load,  compared  well  with  Gu  for  the  material. 

The  results  of  this  analysis  suggests  that: 

1.  The  T*  criteria  is  a  useful  tool  for  evaluating  the  residual  strength  of  damaged 
components. 

2.  Better  experimental  methods  are  needed  in  order  to  determine  more  accurately 
important  orthotropic  laminate  properties  such  as  the  interlaminar  shear  nxxluli 
[63].  Presently,  some  properties,  such  as  Gzj,  arc  difficult  to  determine 
experimentally  and  are  prone  to  experimental  errors. 

3.  More  accurate  and  consistent  methods  of  obtaining  mixed  mode  energy  release 
rates  are  needed.  These  values  would  aid  in  a  more  accurate  and  consistent 
determination  of  failure  loads. 

4.  In  numerical  analyses,  accurate  determination  of  crack  tip  element  stresses  is 
important.  The  use  of  quarter  point  shifts  to  simulate  the  stress  singularity  should 
not  be  ignored,  see  [64]  for  more  details. 

5.  The  effects  of  local  crack  closure  (i.e.,  the  closure  of  delamination  faces)  on  the 
values  of  T*  are  important;  care  must  be  taken  in  finite  element  analyses  to  avoid 
non-physical  solutions  through  the  use  of  appropriate  restraints,  see  [64]  for  more 
details. 

6.  Unlike  cracks  in  metallic  structures  bi-axial  stress  states  may  have  a  very 
significant  effect  on  the  critical  fiacture  (failure)  load.  This  will  be  investigated  in 
more  detail  in  section  7. 


6.  EFFECT  OF  ANISOTROPIC  PLASTICITY 

With  the  development  of  tougher  matrix  materials  such  as  PEEK  the  influence  of 
plasticity  on  the  structural  integrity  of  composite  components  has  become  more 
impxxtant  The  si>read  of  plasticity  and  the  plastic  zone  size  in  the  vicinity  of  the  crack- 
tip  region  are  of  interest  since  it  may  be  adequate  to  employ  fracture  mechanics 
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techniques  if  the  extent  of  plasticity  is  small.  At  present,  data  for  the  plastic  stress-strain 
behaviour  of  composite  materials  are  not  readily  available.  Another  problem  lies  in  the 
definition  of  equivalent  stress  and  strain.  In  this  w<m1c.  the  classical  the(»y  of  anisotropic 
plasticity  is  us^,  together  with  a  reduced  equivalent  stress-strain  curve  experimentally 
determined  by  Chen  and  Sun  [65].  In  this  section  an  incremental  finite-element  analysis 
is  performed  fw  a  three-dimensional  impact  delaminatioa  problem  and  the  effect  of 
plasticity  on  the  fracture  parameter  T*  is  investigated. 

In  the  case  of  incremental  plasticity,  the  value  of  T*  at  the  end  of  a  load  sequence  0. 
P\,...Pn  is  given  by 

N 

r*  =  5^Ar*  (38) 

i»l 

where  AT*  is  given  by 


AT*  =  limf  fAWn, -Afrl^^nds 


and  for  iso-thermal  and  iso-rrtoisture  problems  the  energy  density  W  is 

W  =  (40) 

6.1.  Anisotropic  Plasticity 

This  work  uses  the  classical  three-dimensional  anisotropic  theory  of  plasticity  first 
established  by  Hill.  The  criterion  for  yielding  is  expressed  by  the  'plastic  potential’ 
defined  as 

2/  =  (ar  -O33)'  +^ai.(o33  +|at(oii  +“5(712'  +Y2j'  +Y31') 


•  • 


where 


“t  ~  ^2  ’  “r  “  y2  >  “  jz  ^^^2) 

.K,  T  are  the  tensile  yield  stresses  in  the  fibre  and  transverse  directions,  and  T  is  the  shear 
yield  stress. 

Qearly,  in  the  isotropic  case,  when  X=Y  and  T  =  'Y /-JS,  equation  (41)  reduces  to  the 
familiar  von  Mises  yield  criterion 

2/  =jCtr  [(Ojj  -  0,3 (033  -  a„ )'  +  (o„  -  Oa  )'  -I-  6(7,3"  +  Y23'  +  Y31')]  (43) 

Yielding  is  considered  to  have  occurred  when  the  stress  state  is  such  that  2^1.  In 
classical  incremental  plasticity  theory,  the  plastic  strains  are  defined  as: 

cl£'=^dK  (44) 

do, 

where  dx  is  a  constant  of  proportionality  to  be  determined.  (Carrying  out  the  above 
differentiations,  we  obtain  the  Prandtl-Reuss  equations 
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dEf,  =(0^,0,,  -^0^(0^  -Oj,))d^ 
d£^-iariO^  -O33)- ^<*1.(011 

dE^  =(ar(Oj3  -0a)-^at(0„  -Oa))^  (45) 


dEf,=a.Yu<^ 
dE?,  =a,Yj,dA. 

We  deflne  the  equivalent  stress  a  after  Hill: 

0=  I-,  M.  (46) 

y2(ar+jaJ 

To  evaluate  the  constant  dx  we  define  the  plastic  work  done  as 

dW'  =  ad£^  (47) 


where  is  the  equivalent  plastic  strain  increment.  It  should  be  noted  that  equation  (47) 
is  the  basic  definition  for  ds',  once  a  is  defined.  Another  way  of  expressing  the 
incremental  plastic  work  is 

(48) 

Substituting  equation  (4S)  into  equation  (48)  it  is  ea^y  shown  that 

(M'=ofdK  (49) 

From  equation  (46)  we  see  that 


and  substituting  this  into  equation  (49)  we  have 


Equating  equation  (SI)  with  equaticm  (47)  we  obtain 


Equations  (45),  (46)  and  (52)  completely  define  the  cmhotropic  plasticity  flow  rule  used 
in  this  work.  It  should  be  noted  that  is  estimated  by  Ae  program  from  the  given 
equivalent  stress-strain  curve  and  the  strain  componmts  calculated  using  equation  (45). 


At  yield,  ^1  by  definition,  if  we  substitute  the  values  X=1500MPa  and  T=60MPa  into 
equation  (46)  we  obtain  a  value  for  the  equivalent  yield  stress 

a,  =73.45 MPa  (53) 


4r 
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This  value  of  equivalent  stress  represents  the  limit  of  the  linear  or  elastic  portion  of  the 
stress-strain  curve.  It  should  be  noted  that  the  value  of  X  used  is  the  ultimate  failure 
stress  of  the  composite  in  the  fitne  direction,  since  it  may  be  assumed  that  negligible 
yielding  occurs  in  this  direction  prior  to  failure.  Hence  it  is  clear  that  the  parameter  Oy  is 
significantly  affected  only  by  Y,  the  yield  stress  across  the  fibre  direction.  This 
emphasises  the  assumption  that  plasticity  effects  are  predominantly  dependent  upon  the 
behaviour  of  the  matrix  material. 


Fcmt  the  nonlinear  part  of  the  equivalent  stress-strain  curve,  data  from  wwk  by  Chen  and 
Sun  [65]  are  used.  The  relationship  was  obtained  by  fitting  a  curve  through  experimental 
data.  If  o  is  expressed  in  MPa,  the  relationship  is 


e  =(7.107x10*“)  (o) 


x3.7 


(54) 


Substituting  equation  (S3)  into  equation  (54),  the  equivalent  yield  strain  can  be  obtained; 

Cy  =5.6997x10^  (55) 

An  equivalent  elastic  modulus  may  thus  be  defined; 

(56) 


E  =^=m.S6GPa 

tv 


We  assume  'yielding'  to  have  occurred  when  the  equivalent  stress-strain  behaviour 
surpasses  the  linear  or  elastic  region  of  the  curve.  Another  assumption  used  in  this  work 
is  that  no  account  is  taken  for  the  possibility  that  the  yield  sttess  in  the  across-the-fibre 
direction  may  be  different  in  compression  than  in  tension.  Furthermore,  due  to  the  lack 
of  data  on  shear  yield  stresses,  the  shear  yield  stress  used  here  is  assumed  to  be 

Y 


r= 


V3 


(57) 


although  in  general  this  need  not  necessarily  be  so. 

6.2.  Result  and  Discussion 

The  applied  load  was  increased  in  10  increments  until  the  spread  of  plasticity  was  found 
to  be  extensive.  It  should  be  noted,  however,  that  'plasticity'  for  graphite-epoxy 
composites  is  unlike  plasticity  as  traditionally  understood  for  metals  in  that,  for  the 
former,  there  is  usually  no  cleariy  defined  yield  point.  In  our  simulation,  a  linear  region 
in  the  stress-strain  curve  has  been  artificially  imposed.  Nevertheless,  this  does  not 
invalidate  the  analysis  since  the  object  of  the  analysis  is  to  determine  whether  plasticity 
causes  a  significant  difference  in  the  behaviour  of  T*. 

Figure  14  shows  the  distribution  of  the  ratio  t-c.  the  elastic  value  of  T*, 

and  its  variation  along  the  delamination  front  with  loading  increment  The  maximum 
value  of  T*  is  only  slightly  altered  with  increasing  plasticity.  The  maximum  value  of  the 
ratio  T*pia,tic/T*tia„ie  is  given  in  Table  8.  Significant  plasticity  begins  to  appear  at  load 
increment  4  (40%  of  total  tqrplied  load)  and  it  is  after  this  suge  that  the  deviations  in  T* 
begin  to  become  apparent  (see  Figure  14).  There  are  differences  of  up  to  50%  for  some 
positions  along  the  delamination  front,  although  at  points  of  maximum  T*  (i.e.  at  points 
A  and  C  in  Figure  14)  the  deviations  fiom  the  elastic  solutions  are  very  small.  A  point  of 
particular  interest  is  at  B,  where  with  increasing  plasticity,  T*  changes  sign  from  positive 
to  negative.  In  other  studies,  this  behaviour  has  been  associated  with  local  load-shedding 
[66].  It  should  be  noted  here  that  the  delamination  surface  restraints  applied  to  the  model 


remain  the  same  throughout  the  loading  increments.  These  restraints  are  obtained  from  a 
previous  elastic  run.  To  be  strictly  cc«rect,  however,  the  delaminated  surfaces  should  be 
checked  for  local  crack  closure  for  each  loading  increment  and  the  solution  reiterated.  At 
present,  such  a  procedure  is  unf(»tunately  too  time-consuming  and  costly.  It  is  also 
found  that  the  deviations  of  T*  are  primarily  due  to  the  traction  terms  of  the  integral. 

It  q^>ears  that  the  first  elements  to  become  plasdc  belong  to  the  second  layer  (i.e.  the  45° 
layer)  and  that  plasticity  originates  riom  the  delamination  front  After  load  increment 
number  3  the  plasticity  spreads  in  directions  ±45°  to  the  global  loading  direction.  It 
readily  penetrates  the  0°  layers  above  and  below  the  second  layer  during  and  after  load 
increment  4,  suggesting  that  the  spread  of  plasticity  through  the  thickness  of  composite 
laminates  may  not  be  significantly  restrained  by  the  ply  configuration. 

The  results  of  this  work  indicate  that: 

1.  Plasticity  may  not  significantly  affect  the  severity  of  the  damage,  since  it  does  not 
contribute  to  signiricant  deviations  in  T*^  and  the  strain  energy  density  factor 
[67].  This  implies  that  for  at  least  some  cases,  a  full-scale  plasticity  analysis  is 
not  necessary  and  therefore  signiHcant  cost  in  computing  time  can  be  saved.  For 
the  problem  investigated  here  the  full  plasticity  solution  required  more  than  six 
times  the  CPU  time  of  the  elastic  static  analysis. 

2.  The  7^-integral  applied  to  composite  plasticity  may  become  negative  with 
increasing  plasticity.  It  is  believed  that  this  phenomenon  may  be  caused  by  local 
load-shedding,  see  [66]. 

3.  It  is  shown  that  matrix  plasticity  can  readily  spread  into  adjacent  layers  of 
different  ply  orientations  and  is  not  necessarily  restrained  from  spreading  through 
the  thickness. 


For  metals,  crack-tip  plasticity  can  be  accounted  for  by  replacing  the  crack  length  a  by  an 
effective  length  a+rp  where  rp  is  related  to  the  sire  of  the  plastic  rone.  If  this  approach 
is  adopted  for  the  present  problem,  then  it  is  known  that  once  a  characteristic  damage  size 
has  been  exceeded  the  crack-tip  parameters  are  invariant  to  size  [1],  and  this  would 
suggest  that  the  parameters  should  be  insensitive  to  matrix  ductility.  The  present  work 
appears  to  confirm  this  hypothesis.  The  invariance  of  the  governing  fracture  parameters 
to  plasticity  should  also  hold  for  other  problems  for  which  the  governing  crack-tip 
parameters  asymptote  to  a  constant  value  as  the  size  of  the  crack  increases,  see  [67]  for 
more  details. 


7.  IMPLICATIONS  FOR  FIBRE  COMPOSITE  REPAIRS 

Having  gained  an  undostanding  of  the  mechanisms  governing  the  behaviour  of  impact 
damaged  composites  let  us  now  consider  its  implications  for  repair.  In  this  section 
attention  is  focused  on  the  behaviour  of  impact  damage,  under  both  uniaxial  and  bi-axial 
compressive  loading,  and  a  simple  repair  methodology  is  developed.  This  methodology 
is  then  verified  via  a  simple  coupon  test  program  and  subsequently  by  a  repair  to  an 
impact  damaged  F/A-18  Horizont^  Stabilator. 

There  arc  two  possible  methods  for  repair  impact  damage.  One  approach  is  to  remove 
the  damage  region  and  use  an  internally  bonded  repair.  This  is  very  effective,  but 
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requiics  extensive  bonding  facilities  resulting  in  a  significant  period  of  non-utilization  of 
the  component.  An  alternative  approach  is  to  increase  the  strength  by  reducing  the  net 
sectional  stresses.  This  can  be  achieved  by  placing  an  external  patch  over  the  damaged 
area. 

In  this  section  the  T*  is  applied  to  the  problem  of  an  impact  damaged  composite  structure 
repaired  with  an  externally  bonded  patch.  A  simple  design  methodology  is  presented  and 
substantiated  via  a  laboratmy  test  program  [68].  This  methodology  is  further  illustrated 
by  a  repair  to  a  damaged  htmzontal  stabilator  which  is  tested  to  ultimate  load. 

7. 1.  Numerical  Analysis  of  Bonded  Repairs 

In  order  to  understand  the  mechanisms  governing  failure  and  the  increase  in  residual 
strength  of  a  repaired  structures,  which  contains  impact  damage,  a  detailed  three- 
dimensional  finite  element  analysis  of  the  problem  outlined  in  section  S.l  was 
undertaken.  The  initial  damage  around  the  fastener  hole,  of  9.Smm  diameter,  was 
modelled  as  a  circular  delamination  of  13.75  mm  radius  between  the  second  and  third 
plies  (i.e.  between  the  45°  and  0°  plies).  This  is  an  approximate  simulation  of  the  initial 
damage  where  it  was  found,  from  ultrasonic  C-scanning.  that  the  initial  delamination  was 
nearly  circular,  see  section  5.1.  A  compressive  load  of  150  kN  was  applied  to  the  ends 
of  the  model  in  the  x-direction  (see  Figure  9). 

To  investigate  the  use  of  an  externally  bound  repair  a  12  ply  (0/+45/02/-45/0),  laminate 
bonded  over  the  damaged  area  was  analysed.  This  patch  contained  a  hole  to  allow  for  the 
possible  insertion  of  a  fastener.  To  evaluate  fastener/structure  interaction  effects  various 
boundary  conditions  were  applied  to  the  hole.  These  were  chosen  to  represent:  1)  an 
open  hole.  2)  an  interference  fit  fastener  and  3)  a  bonded  insert. 


7.1.1.  Results  and  Discussion 

For  each  configuration  the  value  of  the  crack  tip  energy  integral  T*  was  computed  at  each 
node  around  the  crack  tip.  A  summary  of  the  maximum  values  of  T*  is  given  in  Table  9. 
The  angles  at  which  the  respective  maxima  occur,  which  arc  referred  to  as  67>  are  also 
shown. 

It  was  observed  that  the  two  plies  above  the  delamination  moved  out  of  plane,  i.e.  crack 
opening  and/or  closure.  Fot  all  of  the  cases  analysed  this  out  of  plane  movement  was 
non-symmetric.  This  asymmetry  results  in  the  asymmetric  growth  of  the  delamination 
which  is  consistent  with  the  C-scans  and  thermal  emission  measurements  [35].  The 
results  of  this  analysis  suggested  that: 

1.  When  the  patch  material  has  the  same  stifihess  as  the  parent  laminate  the 
reductions  in  T*  can  be  estimated  (see  Table  9)  by  multiplying  the  values 
corresponding  to  the  unrepaired  structure  with  the  square  of  the  reduction  in  the 
net  sectional  stress,  see  [69]  for  more  details.  This  infers  that  the  residual  strength 
of  a  repaired  structural  component  can  be  estimated  by  the  following  simple 
formulae: 

Residual  Strength  (repaired)  _  ajunrepaired) 

Residual  Strength  (unrepaired)  o  (repaired) 


•  • 


•  • 
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2.  Prohibiting  local  bending  at  the  hole  had  an  insignificant  effect  on  the  fracture 
parameter.  However,  prohibiting  in-plane  movement  and  local  bending  reduced 
7*.  This  reduction  in  T*  would  ccnrespond  to  an  increase  in  the  failure  load  of 
between  1 1%  to  16%.  As  a  result,  experimental  results,  using  either  a  bonded  or 
an  interference  fit  insert  to  simulate  this  level  of  local  restraint,  are  required  to 
evaluate  this  effect  further. 

3.  Load  bi-axiality  has  a  marked  effect  of  the  fracture  parameter  and  hence  on  the 
failure  of  the  structure. 
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7.2.  Experimental  Program 

72.1.  Specimen  Fabrication 

To  confirm  this  prediction  a  coupon  test  program  was  undertaken.  The  gr^hite  efioxy 
material  used  throughout  these  tests  was  AS4/3501-6  with  a  ply  configuration  of  [(-1-4S2/- 
452A)4)3/90],.  Before  the  specimens  were  cut  from  these  panels,  they  were  C-scanned  to 
determine  the  void  content  of  the  material. 

Each  specimen  was  impacted  with  a  12mm  (l/2“)  diameter  ball  bearing  with  a  mass  of  1 
kg  and  from  a  height  of  1.3m.  A  summary  of  the  absorbed  energies  is  listed  in  Table  10. 
The  extent  of  damage  was  limited  by  the  use  of  20,  30  and  40mm  "windows"  clamped 
around  the  specimen  during  the  impacting.  For  the  single  impact  specimens  all  three 
diameter  "windows"  were  used  to  produce  damage  with  diameters  ranging  from 
approximately  19.8tnm  to  39.7mm.  A  6mm  diameter  hole  was  also  drilled  through  the 
centre  of  the  damaged  region  of  these  specimens.  All  specimens  were  subjected  to  a 
C-scan  of  the  impacted  area  and  the  damage  size  was  approximated  from  these  scans.  A 
trxne  detailed  description  is  given  in  reference  [68]. 

7.22.  Repair  Fabrication 

In  order  to  validate  the  simple  design  rule  previously  postulated,  it  was  required  that  the 
effective  stiffness  and  ply  configuration  of  the  patch  be  representative  of  the  parent 
material.  The  material  us^  fOT  the  patch  was  AS4/3S01-6  and  was  16  plies  thick,  with  a 
ply  configuration  of  [02/+45/-45/-45/+45A)2],.  The  ends  of  the  patch  were  scarfed  to 
reduce  the  peel  stresses  in  the  adhesive,  see  [^].  The  length  of  the  patch  and  distance  to 
the  edge  of  the  grips  were  190mm  and  60mm,  respectively.  The  patches  were  bonded  to 
the  parent  laminate  using  the  cold  setting  acrylic  adhesive  FLEXON  241.  This  adhesive 
was  chosen  for  its  excellent  shear  strength,  ease  of  application  and  because  environmental 
effects  were  not  an  issue  in  diis  lest  series. 

723.  Test  Methodology 

Unrepaired  specimens  were  tested  individually  while  repaired  specimens  were  tested 
back-to-back  [68].  Four  strain  gauges  were  bonded  on  the  unrepaired  specimen;  two 
above  and  two  below  the  damage.  Each  pair  was  positioned  on  the  mid-width  and  on 
opposite  faces  in  order  to  determine  axial  and  bending  strains.  The  repaired  specimen 
only  had  two  strain  gauges  located  on  the  patched  side.  In  each  case  the  gauges  were 
1 10mm  from  the  centre  of  the  specimen.  The  specimens  were  loaded  in  compression  and 
to  avoid  the  problem  of  global  buckling  an  anti-buckling  rig  was  used  [68]. 


•  • 


•  •  • 
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7.2.4.  Results  and  Discussion 

All  specimens,  except  two  which  exhibited  extensive  bending,  produced  load  versus 
strain  curves  which  were  essentially  linear  to  failure.  Specimens  4  and  5  exhibited 
extensive  bending  near,  and  up,  to  failure  thereby  reducing  the  failure  load  obtained.  The 
failure  strains  for  each  specimen  are  tabulated  in  Table  10.  These  results  show  that  the 
failure  strains  follow  the  asymptotic  nature  outlined  in  section  2.  The  variation  in  the 
failure  strains,  f<x  a  given  damage  size,  reflects  variations  in  specimen  geometry  as  well 
as  variations  in  the  structure  of  the  internal  damage  due  to  impact.  In  all  cases,  the 
patches  and  the  adhesive  bond  failed  after  the.  parent  laminate. 

The  failure  load  of  the  repaired  specimen  can  also  be  predicted  (see  Table  10)  using 
equation  (58)  and  only  requires  a  knowledge  of  the  unrepaired  residual  strength  and  the 
change  in  net  sectional  stresses  due  to  patching.  The  change  in  the  net  sectional  stress 
can  be  calculated  directly  from  the  change  in  net  sectional  area.  This  result  is  believed  to 
substantiate  the  trends  predicted  in  the  previous  .section  and  significantly  simplifies  the 
repair  design  philosophy. 


8.  DAM  AG  ED  ASSESSMENT  OF  AN  F/A- 18  HORIZONTAL  STABILATOR 

8.1.  Background 

As  an  example  of  how  this  approach  can  be  used  to  repair  damaged  suuctural 
components  let  us  consider  the  repair  of  a  damaged  F/A- 18  horizontal  stabilator,  which 
has  recendy  been  incorporated  into  the  (Domposite  Repair  Engineering  Development 
Program  (CTREDP).  CREDP  is  a  joint  program  between  the  Canadian  Forces  (CF),  the 
RAAF  and  the  United  States  Navy  (USN)  to  evaluate  the  repair  capability  of  damaged 
composite  components  on  the  F/A- 18.  As  part  of  the  CREDP  program,  the  Aeronautical 
Research  Laboratory  (ARL)  was  tasked  to  develop  and  analyse  a  repair  for  a  damaged 
horizontal  stabilator,  which  was  initially  classed  as  unserviceable  and  irreparable  due  to 
two  fragment  strikes  from  a  tracer  rocket. 

The  F/A- 18  horizontal  stabilator  comprises  an  aluminium  honeycomb  sandwich  structure 
with  graphite/epoxy  skins,  with  fibres  cniented  in  the  0°  and  ±45  directions.  The  cross- 
section  of  the  stabilator  is  fully  symmetric.  In  the  region  of  the  damage  the  skin  is  29 
plies  thick  (3.68mm),  but  tapers  off  to  7  plies  (0.89nim)  at  the  leading  edge  and  has 
incurred  extensive  local  damage  to  the  composite  sidn  and  underlying  honeycomb,  but 
neither  fragment  penetrated  to  the  other  side  of  the  stabilator.  As  shown  in  Figure  15,  the 
two  damaged  areas  are  approximately  the  same  size;  however  there  was  a  slight 
difference  in  the  impact  angle  of  the  two  fragment  strikes.  Hie  location  of  the  danutged 
zone  and  ply  orientation  can  be  seen  in  Figure  16. 

8.2.  Structural  Evaluation 

In  order  to  evaluate  the  structural  signiflcance  of  the  damage,  i.e.  if  the  original  design 
was  optimised  for  impact  damage,  a  test  rig  capable  of  ^plying  the  design  loads  was 
used.  The  stabilator  was  mounted  in  the  rig  by  means  of  the  spindle,  which  is  used  to 
connect  the  stabilatOT  to  the  aircraft.  Hence  all  bending  loads  applied  to  the  stabilator  are 
transferred  to  the  rig  via  the  spindle.  A  lever  arm  connected  to  the  root  of  the  stabilator 
transfers  all  the  torque  loads  to  the  test  rig.  Static  loads  were  applied  to  the  structure  via 
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two  flexible  air  bags.  The  design  of  the  test  rig  is  such  that  it  allows  these  air  bags  to  be 
mounted  above  or  below  the  test  article,  as  described  in  reference  [70]. 

Four  strain  gauge  rosettes  were  located  in  a  rectangular  pattern  around  the  damage  zone 
and  in  the  ccHiesponding  location  on  the  undamaged  sldn,  see  Figure  16.  The  gauges  in 
each  rosette  were  aligned  in  the  direction  of  the  0  ami  ±45°  fibres.  Two  displacement 
transducers  were  placed  halfway  between  the  left  hand  inboard  rosettes  (1  and  2),  in  the 
0°  fibre  direction,  as  seen  in  Figure  16.  The  gauge  length  of  the  displacement  transducers 
was  38Smm.  The  strain  gauge  rosettes  and  di^lacement  transducers  were  positioned  so 
as  to  evaluate  the  change  in  compliance  of  the  structure  due  to  the  damage. 

Displacement  and  pressure  transducers  were  required  to  measure  the  extension  and 
pressure  of  each  air  bag.  The  force  applied  by  each  air  bag  is  a  function  of  these  two 
variables  and  was  determined  from  a  calibration  graph.  The  root  bending  moment 
(RBM)  was  calculated  and  used  as  the  reference  load  applied  to  the  structure.  Three 
extra  displacement  transducers  were  attached  to  the  stabilator  at  the  tip,  leading  and 
trailing  edges,  to  measure  tip  displacement  and  torque  induced  by  the  applied  load. 

82.1.  Test  Description 

The  first  stage  of  the  investigation  involved  the  application  of  a  static  load  by 
incrementing  the  air  bag  pressure  in  steps  of  10  kPa.  Here  the  maximum  load  applied  to 
the  structure  was  not  considered  to  be  critical  since  the  primary  objective  was  a 
comparison  between  the  strains  and  the  compliance  on  the  top  and  bottom  surfaces  of  the 
stabilator.  Initially,  the  air  bags  were  placed  underneath  the  stabilator,  producing  tension 
in  the  damaged  surface  skin,  and  a  dummy  loading  and  unloading  run  was  performed  to 
allow  the  structure  to  settle.  Several  loading  artd  unloading  runs  were  made  in  this 
configuration  with  strain  gauge  and  transducer  readings  being  conducted  at  each 
increment  of  pressure.  The  loading  was  then  repeated  with  the  air  bags  re-configured 
above  the  stabilator  producing  a  compressive  load  in  the  damaged  surface. 

822.  Results  and  Discussion 

The  Uldmate  Bending  Moment  G-^BM)  and  the  Design  Limit  Bending  Moment  (DLBM) 
for  the  horizontal  stabilator  ate  120.3  kNm  (1065  kip-in)  and  80.2  kNm  (710  kip-in), 
respectively  [70].  The  RBM  was  calculated  for  each  load  step  and  the  test  achieved  31% 
(47%)  and  26%  (39%)  of  the  UBM  (DLBM)  in  tension  and  compression,  respectively. 

The  strain  survey  results  for  one  loading  cycle  are  shown  in  Table  1 1,  for  the  case  when  a 
compressive  load  was  applied  to  the  damaged  surface.  The  results  for  the  case  of  tensile 
loads  iqiplied  to  the  damaged  surface  are  very  similar.  All  strain  values  are  presented  in 
miCTOstrain.  Table  1 1  shows  that  there  was  no  structurally  significant  difference  between 
the  strain  gauge  results  for  strains  in  the  0°  flbre  direction  on  the  top  (damaged)  and  the 
bottom  (undamaged)  surfaces,  except  for  gauge  D2  which  has  a  24%  to  27%  reduction 
depending  upon  the  value  of  the  RBM,  due  to  damage.  This  may  be  due  to  the  gauge 
being  shielded  by  the  damage.  The  compliance  readings,  tabulat^  at  the  bottom  of  the 
tables  as  U'  for  Ae  undamaged  surface  strain  and  *D'  for  the  damage  surface  strain,  show 
no  significant  difference  between  the  top  and  bottom  surfaces.  This  implies  that  the 
global  compliance  has  not  been  affected  by  the  damage  and  indicates  that  there  is  little 
structural  degradation.  This  was  confirmed  by  a  thermal  elastic  evaluation  of  the 
damaged  region,  see  [71]. 
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8^  The  Finite  Element  Model 


To  understand  the  mechanisms  involved  in  this  problem  and  to  assist  in  designing  a 
repair  a  finite  element  analysis  was  undertaken.  The  damage  section  of  the  horizontal 
stabilatOT  was  modelled  by  representing  the  skin  as  a  two-dimensional  membrane  with 
effective  laminate  properties  and  the  underlying  honeycomb  as  three-dimensional  iso¬ 
parametric  iHick  elements.  The  small  changes  in  the  compliance  and  the  local  strains  due 
to  damage,  as  observed  in  the  previous  section,  suggest  the  possibility  of  using  an 
external  doubler  to  repair  the  damaged  stabilator.  In  this  case  the  adhesive  layer  was 
modelled  using  three-dimensional  iso-parametric  tmck  elements  with  the  doubler 
represented  by  two-dimensional  memlnane  elements,  with  an  effective  laminate  property. 
Tluee  cases  were  analysed,  the  undamaged  sldn,  the  damaged  but  unrepaired  skin  and  the 
damaged  and  repaired  skin.  The  damaged,  but  unrepaired,  model  contained  520  nodes, 
128  elements  and  1255  degrees  of  freedom  and  the  repaired  nxxlel  contained  818  nodes, 
256  elements  and  1976  degrees  of  freedom. 

The  material  used  for  the  skin  and  the  doubler  was  again  AS4/3501-6  and  the  associated 
laminate  properties  can  be  found  in  Table  12.  Young's  modulus  of  the  honeycomb  and 
the  adhesive  was  taken  to  be  2.0  GPa  and  1.8  GPa,  respectively,  with  a  Poisson's  ratio  of 
0.35  for  both.  The  doubler  was  16  plies  (2.032  mm)  thick  and  a  pressure  of  400  MPa 
was  applied  along  one  edge  of  the  skin. 

The  accuracy  of  the  finite  element  model  was  evaluated  by  comparing  the  predicted  and 
measured  strains,  obtained  in  [70,72].  For  the  undamaged  stabilator  the  ultimate  design 
bending  moment  corresponds  to  a  design  limit  strain  of  approximately  4850  |i£,  see 
[70,72].  The  corresponding  predicted  strains  at  the  two  trailing  edge  strain  gauges,  1  and 
3,  were  approximately  4950  pe  and  4860  pe,  respectively,  and  the  prediaed  value  of 
strain  over  the  region  measured  by  the  LVDT  was  also  close  to  the  measured  value. 
However,  because  the  actual  structure  tapers  off  towards  the  leading  edge,  the  strains  in 
the  numerical  model  closest  to  the  leading  edge  (gauges  2  and  4)  were  higher  than  those 
achieved  in  the  experiment  Nevertheless  it  was  felt  that  this  model  would  be  useful  in 
predicting  the  change  in  strain  due  to  the  application  of  an  external  doubler. 

83.1.  Predated  Reduction  in  Strain  Due  to  on  External  Doubler 

Following  the  development  of  the  undamaged  finite  element  model,  two  cases  were 
considered.  The  first  case  simulated  the  damaged,  but  unrepaired,  stabilator  and  the 
second  case  modelled  a  bonded  external  doubler.  A  summary  of  the  predicted  strains 
measured  by  the  LVDT  and  the  strains,  in  the  0°  fibre  direction,  at  gauges  1-4  are  given 
in  Table  13. 

The  application  of  the  external  doubler  significantly  reduced  the  strain  around  the 
damaged  zone  leading  to  an  average  reduction  in  strain  of  28.5%.  The  peak  adhesive 
shear  stress  was  less  than  25  MPa  and  was  below  the  endurance  level  of  the  thin  film 
adhesive  FM73,  see  [73],  indicating  that  the  repair  should  not  debond  and  should  be 
inherently  damage  tolerant,  see  [73].  The  stress  concentration  around  the  damage  was 
also  found  to  be  very  localised  and  was  consistent  with  the  experimentally  measured 
thermal  emission  profile  [72],  as  well  as  with  the  experimentally  measured  values  of 
strain  and  compliance. 


8.4.  Repair  Evaluation 


Fc^wing  the  initial  tests  on  the  htxizontal  stabilator  the  damaged  region  was  repaired 
using  an  external  doubler  which  was  tailored  to  matched  the  stifhiess  of  the  stabilator 
sldn.  The  ply  configuration  of  the  repaired  laminate  was  (+45/-45A)2/9Q/-45/+45A)),. 
The  loose  fibres  around  the  impact  zone  were  removed  and  the  edges  of  the  damage  zone 
were  scarfed.  The  whole  region  was  then  filled  with  pacldng  adhesive  and  the  repair 
laminate  was  then  bonded  onto  this  surface  using  FM73  adhesive. 

The  test  procedure,  outlined  in  section  8.2.1,  was  then  repeated  with  additional  strain 
gauge  rosettes  placed  on  the  repair.  The  location  of  the  repair  ranged  from  0.610m  (24") 
to  1.092m  (43")  from  the  root  of  the  stabilator.  The  repair  successfully  withstood  both 
tension  and  compression  loading,  see  [74].  In  tension,  the  majority  of  the  region 
exceeded  the  DLBM.  The  measured  strain  reduction  was  approximately  28%  and  was 
equivalent  to  the  reduction  in  the  nett  section  stress. 

8.5.  Remarks 

This  section  has  presented  a  repair  methodology  that  can  be  used  as  a  first  approximation 
for  field  repairs  to  damaged  composite  structures.  It  has  also  been  shown  that  an 
externally  bonded  doubler  was  able  to  repair  impact  damage,  even  for  thick  structural 
components.  The  use  of  externally  bonded  patches  has  the  advantage  of  being  quick  and 
easy  to  apply.  The  F/A-18  horizontal  stabilator  repair  also  indicated  that  extensive 
removal  of  damaged  material,  excluding  penetration  clean  up,  can  be  avoided  thereby 
making  the  repair  process  particularly  simple. 


9.  CONCLUSIONS  AND  RECOMMENDATIONS 

This  report  has  attempted  to  clarify  the  fundamental  mechanisms  responsible  for  the 
degradation,  in  both  sutic  strength  and  fatigue  performance,  due  to  low  velocity  impact 
damage.  It  has  been  shown  that  for  the  development  of  a  damage  tolerant  philosophy  a 
consistrat  test  and  analysis  methodology  is  required  for  measuring  and  computing  the 
critical  energy  release  rates.  For  Modes  11  and  lU  this  methodology  should  allow  for  both 
in-elastic  effects  and  for  strain  rate  effects. 

It  has  also  been  shown  that  for  delamination  damage  in: 

(a)  composite  laminates, 

(b)  at  step  lap  joints,  or 

(c)  at  mechanically  fastened  composite  joints, 

a  stage  is  reached  after  which  a  significant  increase  in  the  size  of  the  damage  does  not 
significantly  reduce  the  residual  strength.  This  significantly  simplifies  the  methodology 
for  estimating  critical  damage  size.  It  has  also  been  shown  that,  in  fatigue,  there  is  a  well 
defined  endurance  limit  and  that  small  changes  in  the  cyclic  stress  amplitude  can  produce 
very  large  changes  in  fatigue  life.  In  contrast  to  the  situation  for  cracks  in  metallic 
structural  components,  load  bi-axiality  has  a  marked  effect  on  residual  strength.  This 
means  that  for  sovice  aircraft  the  residual  strength  and  the  strength  after  repair  cannot  be 
estimated  using  either  uni-axial  test  specimens  or  via  simplistic  one  dimensional  analysis. 
This  conclusion  needs  to  be  validated  experimentally. 
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A  simple  repair  design  methodology,  that  can  be  used  as  a  first  ^proximation  for  field 
repairs  to  damaged  composite  structures,  has  also  been  presented  and  it  has  been  shown 
thiu  an  externally  bonc^  doubler  was  able  to  repair  impact  damage,  even  for  thick 
structural  con^xxients.  The  F/A-18  Htaizontal  Stabilattx'  repair  indicated  that  extensive 
removal  of  damaged  material,  excluding  penetration  clean  up,  can  be  avoided  thereby 
making  the  repair  process  particularly  simple. 
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TaMe  la:  CompariaoB  of  experimental  values  of  Gi;  from  the  literature 


Mattrial 

Test  method 

Gic 

(JInP) 

T300/5208 

Wilkins  efai.'* 

DCB 

87.6 

T300/5206 

Ramkuinar  &  Whitcomb’’ 

DCB 

102.6 

1300/5208 

Ouu’^ 

DCB 

86±8 

T300/5208 

CrBrien23 

EDT 

137 

T300/5208 

Byers“ 

Modified  width-tapered  DCB 

160-241 

GY70/5208 

Russel  &  Street^o 

DCB 

83 

ASl/3501-6 

Russel  ic  Street*’-^ 

DCB 

no 

ASl/3501-6 

Jurf  &  Pipes^ 

Modified  Arcan  Test 

79 

AS4/3501-6 

Aliyu  &  Daniel^’ 

DCB 

198-232 

ASl/3502 

Bradley  &  Cohen’^ 

DCB 

155 

ASl/3502 

Whitney  &  Browning** 

DCB 

140 

ASl/3502 

NkhoUs  6c  Gallagher^ 

DCB 

105-175 

ASl/3502 

Bradley  6c  Cbhen’^ 

CT  (Omipact  Tension) 

225 

ASl/3502 

Whitney  6c  Browning** 

90“  CN 

154 

ASl/3502 

Whitney  6c  Browning** 

EDT 

267 

AS4/3SQ2 

Bradley  6c  Cohen’* 

DCB 

225* 

AS4/3502 

Bradley  &  Grfien’* 

CT 

120 

AS4/3502 

Whitney  &  Knight*^ 

EDT  SOB  mm  wide  sptecimen 

103 

AS4/3502 

Whitney  6c  Knight** 

EDT  25.4  mm  wide  specimen 

238 

AS4/3502 

Whitney  6c  Knight** 

EDT  modified  with  inserts 

81-% 

AS4/3502 

Whitney  6c  Knight** 

DCB 

158 

T300/1034C 

Donaldson** 

90“  CN 

77.9 

T300/934 

O'Brien** 

EDT 

216 

T300/934 

Wang  et  al.*’ 

DEN  (Double  Edge  Notched) 

157 
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Tabte  lb:  Description  of  analysis  and  experiment  for  determining  Gfc 

^tference  Teat  Description 

metitod 

Wilkins  ef  a/.**  DCB  Compliance  method  of  test  procedure  em^rfoyed,  i.e.  G  = 

P^(dClda)l2W.  Load-dispUcement  curve  nv>nitofed  for  static 
loading.  Relationship  between  orntpliance  and  crack  length  was 
obtained,  C  =  2a^/3El  =  Also,  the  relationship.  Pc  = 

(GcWED^I^Ia  =  Alia  was  obtained.  The  value  of  Gc  can  then  be 
calculated.  Both  experimental  curves  of  C  vs  a  and  Pc  vs  a  were 
obtained  via  a  least-squares-fit  routine  to  obtain  A]  and  Aj. 

Ramkumar  &  DCB  Similar  data  analysis  as  in  Wilkins  et 

Whitcomb'® 


Russell  &  Street^  DCB  Similar  data  analysis  as  in  Wilkins  et  alJ^. 


Bradley  &  DCB 
Cohen'^ 


Aliyu  &  Daniel^'  DCB 


Whitney*  DCB 
Browning^ 


Chai'^  DCB 


Analysis  based  on  linear  beam  theory.  G/  =  8Pj/L^/WEI,  where  P,  is 

the  applied  load,  te  is  the  length  of  the  split  laminate.  The  analysis 
was  essentially  similar  to  Wilkins  et  alJ^. 

Beam  analysis  method  used.  The  analysis  was  extended  to  include 
transverse  shear  deformation  and  kinetic  energy  effects  on  the 
energy  release  rate.  A  semi-empirical  method  based  cm  tiie 
assumption  of  partial  elastic  foundation  support  beyond  the  crack 
tip  was  propos^.  Effect  of  crack  velocity  on  energy  release  rate 
studied.  Gfc  found  to  increase  with  crack  velocity. 

Gfc  obtained  duough  the  relationship  Gfc  =  l/26Ac(/*,82-P25,). 
where  Pj  and  P2  are  the  loads  before  and  after  crack  extension,  8j 
and  82  the  deflections  before  and  after  crack  extension,  and  Ac  the 
increment  in  crack  extension.  The  analysis  was  based  on  area 
between  the  loading  and  unloading  curves. 

Semi-empirical  method  of  analysis  used.  An  empirical  relationship 
between  compliance  C  and  crack  length  1,  C  =  a/",  where  a  and  n  are 
constants,  was  used  to  obtain: 


Jurf*  Pipes?®  Modified  The  aluminium  fixtures  enabled  mixed-mode  loading  of  the 
Arcan  test  specimen  by  varying  the  angle  of  loading  a.  For  pure  Mode  I 
^”8*®  loading,  a  =  90”  and  the  an’Ued  stress  normal  to  the  crack  plan  was 
<^  =  f*/A,  where  Pis  the  applied  load  and  A  is  the  area  of  qredmen 
Notch  on  which  the  load  is  applied.  There  was  some  scatter  in  the  values 
of  Xfc  obtained.  Gfc  was  obtained  fiom  Xfc. 
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TaUc  lb:  Coot'd 


RtfmHce 

Test 

method 

Description 

Byw“ 

DCB 

Width-tapered  speciniens  fm  comparison  of  interlaminar  fracture 
toughness  of  composites  with  diffemt  resin  systems.  Ho%vever,  dw 
aluminum  tabs  were  actually  half-inch-th^  l^tes  adhesively 
bonded  to  the  entire  top  and  bottom  surfaces  of  the  lamiiuited 
spedmen.  The  values  of  Gfc  obtained  were  much  higher  than  diose 
obtained  by  others  using  normal  DCB  specimens,  probably  due  to 
increased  stiffness  introduced  by  the  aluminium  plates.  The 
expression  Gfc  =  12J’^/fEAX«/bl^  was  used,  where  P  is  the  crack 
initiation  load  and  (o/b)  is  die  specimen  taper  ratio. 

Donaldson^ 

90“  CN 

Gfc  was  obtained  through  calculation  of  =  Offc-Jiu,  where  a  is 

the  half  crack  length,  and  o^c  ^  far-field  stress  normal  to  d>e 

crack.  G/  was  associated  with  X/  via  an  equation  containing  d\e 
orthotropic  material  properties. 

Whitney  &: 
Browning^ 

90“  CN 

Gfc  obtained  through  calculation  of  Ki^  =  a“^n(a+aQ),  where  Og  is 
an  inherent  flaw  determined  from  the  average  stress  criterion,  o~  is 
the  notch  strength  of  plate  of  infinite  extent.  was  obtained  from 

the  strength  of  a  tensile  coupon  with  dnite  width,  via  an 
isotropic  width  correction  relationship.  Gfc  was  obtained  from  Xfc. 
Gfc  obtained  was  higher  than  values  obtained  dirough  0“  DCB  tests. 

Wang«tai3i 

DEN 

G(  available  at  notch-tip  calculated  by  a  two^limensional  plane- 
stress  finite  element  routine. 

Bradley  & 
Cohen’* 

CT 

Displacement-load  relationship  monitored.  Tests  conducted  in 
displacement  control  to  allow  stable  crack  growth.  Standard  linear 
fracture  mechanics  relationship  G/=(P^/2B)(dC/da)  used. 

aBrien23-2* 

EDT 

Laminated  plate  theory  and  simple  rule  of  mixture  for  stiffness  loss 
due  to  delaminations  were  used  to  derive 

zIt 

where  Elam  is  stiffness  of  plate  calculated  from  laminated  plate 
theory,  £*  is  the  stiffness  of  pnlially  delaminated  laminate,  is  the 
critical  strain  at  the  onset  (X  delamination,  and  T  is  the  laminate 
thickness. 

4r 


•  • 


•  • 


•5/ 


Table  lb:  Coat'd 
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Refermce  Test 
method 


Description 


Modified  Specimens  had  hee-edge  starter  cracks  of  length  a  embedded  in  a»e 
Knight"  EDT  centre  plane  of  specimens.  Classical  laminated  plate  thewy  invoked 
to  obtain  G,,  =  Ae*(£, -£•).  where  k  is  the  half  duckneas  of 
laminate,  is  the  critical  strain,  and  £.  and  £*  are  modified  inplane 
stiffness  terms.  E*  here  is  not  the  same  as  O’Brien's  £*.  Values  of  Gfc 
obtained  slightly  lower  than  those  from  unnvodified  EDTs. 


Whitney  & 
Browning^ 


EDT  Laminated  plate  theory  used  to  derive  =te2(l-£,*  /£,)£i  where 
h  is  half  latrrinated  thickness,  is  the  critical  strain  at  delantirration, 
Ej  is  the  laminate  modulus  in  the  load  direction,  £,*  is  an  effective 
laminate  modulus  for  delaminated  composite,  aiui  £i  is  the 
experimentally  determined  value  of  Ej.  Both  E}  and  £*  were 
calculated  from  laminated  plate  theray.  If  £i  =  Ej,  Aen  Ae 
expression  for  Gfc  would  be  the  same  as  O'Brien's  expression. 
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TaMe  2:  Comparison  of  experimealal  values  of  Guc  from  ihe  iiieiatiHe 

Rtfemtce  Description  of  test  method  and  analysis  Material  Cuc 

Wilkins et  al.^^  CLS  (Cracked  Lap  Shear)  specimens:  Relative  proportions  T330/5206  154 

of  Gf  and  Gn  (Stained  through  riniie  elemmit  aitalysis. 

Proportion  of  Mode  II,  G  found  to  be  approx.  75%  cf  total 
energy  release  rate. 

Rainkuinar&  CLS  specimens.  The  embedded  delamination  did  irot  T300/5206  876,criterion(a) 
Whitcombe^’  always  propagate  along  the  midplane,  but  crossed  over  to  456,criterion(b) 

the  adjacent  interface.  A  simple  strengtir-of-materials  643,criterion(c) 

analysis  of  CLS  specimen  to  account  for  the  change  in 
stiffiress  due  to  tl^  crossover  was  used  to  compute  the 
change  in  compliance  with  crack  length,  dC/da.  A 
geometrically  iron-luwar  finite  element  analysis  was 
performed  to  determiite  the  relative  proportions  of  Gi  aivl 
G//  at  the  critical  load.  Gc  was  obtained  experimentally. 

Values  of  G;,  G/;  and  Gfc  were  substituted  into  three 
failure  criteria  and  G(ic  obtairred: 

(a)  G//Gfc+G;//G/fc=l, 

(b)  (G,IGuP^G„IG,k)^=h 

(c)  (G,IGk^+(G„IG,kP*(GilGk)+(Gn/G,ic)=l. 

Russell  &  ENF  (EtKl  Notched  Flexural)  specinwns.  Specimerrs  were  ASl/3501-6  452+6 

Street^®'^^  subjected  to  three  point  flexural  loading.  The  Teflon  HMS/3501-5  152+11 

started  rrotch  tip  was  placed  midway  between  two  AS4/2220-3  75()+25 

loading  pins.  Beam  theory  aruilysis  was  used  to  obtain 
dCjda  to  be  substituted  into  the  expression  for  G/fc, 
yielding  Gitc=9a^P^Cbf2b(2L^+3(^),  where  Cs  is  dhe 
compliance  due  to  bending. 


Chatterjeert  Three  Point  Bend  Tcjt  specimens.  Implartted  defects  in  ASl/3501-6  1042+214 

the  form  of  two-ply  Teflon  delaminatiorts  were  located  in  Blunt  crack 

the  region  of  compressive  flexural  stress.  A  computer  tips 

code  was  used  to  ^ve  the  problem  of  disbonds  between  95(>+175 

two  piles.  Some  of  the  specimens  were  fatigued  to  Sharp  crack 

produce  sharp  crackfronts.  tips 


Donaldson^^  Modified  Three  Rail  Shear  Test.  This  was  rrot  a  T300-1034C  506 

delamination  test,  but  an  in-platre  crack  was  cut  into  an 
urudirectiorral  laminate  parallel  to  the  flbre  direction. 

Earlier  flnite  element  arulysis  by  Lakshmirrarayana 
showed  a  uruform  state  of  in-plane  shear  stress 
throughout  the  central  region  of  the  specimen.  This  shear 
stress  was  computed  via  ic=Pf/2il,  where  Pc  is  the 
critical  applied  load,  t  and  I  are  the  specimen  thickness 
and  len^  respectively.  The  stress  intensity  factor 
A/fc  =  TcViM  where  a  is  the  half  crack  lengtii,  was 
computed  rnd  die  value  of  Goc  thus  calciilated. 


-36- 


Table  2:  Cont'd 

Referrnce  Description  of  test  method  and  analysis  Material  Gik  & 

(ihnh 

JurffcPipes^  Modified  Arcan  Test  on  Single  Edge  Notch  specimens-  ASl/3501-6  670 

The  fixtures  enabled  mixed-mode  loading  of  the  specimen 
by  varying  the  angle  of  loading,  a.  The  stresses  applied  to 

the  specimen  were  simply  a  =  (Psiaa)/A  for  the  normal  • 

stress,  and  t  =  (f*oosa)/A  for  the  shear  stress,  where  A 

was  the  area  of  specimen  through  which  the  load  was 

transmitted.  The  specimens  were  91  ply  thick,  with  teflon 

inserts,  and  where  adhesively  bond^  to  the  aluminium 

Arcan  fixtures.  For  Mode  II  loading,  a  =  0°.  Some 

specimens  failed  at  the  metal-composite  adhesive  bond-  * 

K/k  was  obtained  and  related  to  Gnc- 


Table  3:  Summafy  of  experimeaial  data  oa  S-N  curves  derived  bom  the  literature 


R^mtux 

Maximum 
eompresewe 
stnin  range 

Maximum 
compressive 
stress  range 

Range  of 
cycUs  to 
failure 

Comments 

Rosenfeld  & 
Gavtae'^ 

2400-3000|t 
(i.e.  de  z 
600|ie) 

2x10^-2x10^ 

Laminated  plate  spedmens; 
AS/3501-6;48ply; 
Delantination  area  12.9  cm^; 
16.4  J  impact;  R  =  -<*‘ 

Demuts  et  al*^ 

Strain  fixed  at 
3000|ie 

10* -10^ 

Multirib  specintens; 
AS6/2220-3; 

136  J  impact;  R  =  10. 

Demuts& 

Horton^ 

S  =  0.6^.7* 

5x10* -10^ 

3-stiffener  panel; 

AS6/2220-3; 

136  J  impact;  R  =  10. 

Walter  ef  al.*^ 

S  =  0.28-0.30 

10* -10^ 

Laminated  plate  specimens; 
T300/5208; 

6.2  J  impact;  R  =  -1.0. 

Stellbrink  & 
Aoki^* 
Stellbrink** 

200-250  MPa 
(i.e.  Aa  = 

50  MPa) 

0-105 

Laminated  plate  specimens; 
T300/CODE69;  18  ply; 
Impacted  with  ball  of  mass 
255  g  and  radius  of  5  mm  at 
velocity  of  6.15  m/s;  R  =  -1.0. 

Potter^ 

Stress  level 
fixed  at  206 
MPa 

622  -  8242* 

Tapered'thidcness  spedmer«; 
XAS/914;32ply-24ply; 

5.0  J  impact;  R  =  10. 

Ryder  et  ai.^ 

Stress  level 
fixed  at  241 
MPa 

Stress  level 
fixed  at  152 
MPa 

2x10* -lO^ 

2x10*-5x105 

24  ply  damaged  hole 
specimens; 

R  =  -1.0. 

32  ply  damaged  hole 
specimens; 

R  =  -1.0. 

Sis  the  ratio  of  the  average  peak  foUgue  load  to  damaged  residual  static  strength. 
Test  terminated  before  spedmen  failed. 
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Tabit  4:  Values  of  ihe  exponent  V  in  the  growth  law  (daldN)»c(Gf  as  predicied  by  venous 
studies  for  Mode  I,  Mode  II,  and  Mixed  Mode  cases 


Reference 

Model 

Exponent  ‘n' 

Mixed-Mode  or  Mode  II 

Wilkins  e(  a/.'” 

Wang  et  al.*^ 

20-30 

7-8  (Mixed-Mode) 

Rainkuinar  &  Whitcomb^ 

8-10 

6  (Mixed-Mode) 

O’Brien" 

3-3J 

. 

Chatteijee  et  al.*^ 

- 

2.4  (Mode  0) 

Whitcomb*® 

15-20 

- 

Table  5:  Values  of  the  exponent  V  in  the  growth  law  <da/d/V)=cf  A  G;//*  as  developed  by 
Russell  h  Street^^  in  foe  order  of  increasing  toughness 


Exponent  *n' 

Material 

5.79 

ASl/3501-6 

First  generation  high  temperature  graphite-epoxy 

5.71 

AS4/2220-3 

Second  generation  high  temperature  graphite-epoxy 

4.52 

C6000/n55 

An  intermediate  temperature  rubber  toughened  graphite-epoxy 

3.88 

AS4/APC2 

High  temperature  graphite-thermoplastic  composite 

Table  <S:  Properties  of  T30Q/S208 


Ell 

E22 

VJ2 

Gi2 

Gj3 

G23 

(GPa) 

(GPa) 

(GPa) 

(GPa) 

(GPa) 

• 

• 

Case  1  -  Equal 
shear  moduli 

141.33 

945 

0.31 

5.18 

5.18 

5.18 

Case  2  -  Unequal 
shear  moduli 

141J3 

945 

0.31 

5.18 

5.18 

3.63 

• 

• 

TABLE  7:  Properties  of  nuitrix  nuiterial 


Em  'Om 

(GPa) 


350 


0.3 


TabkS 


Uamtmm  value  <4  the  ratio  plastic  to  elastic  fracture 
parameter 


^^=L157 


Tabic  9:  Summary  of  T*  (max  value)  and  its  angular  location  6^ 


Case 

T* 

0  T 

Predicted  T* 

lUN 

0.0907 

30* 

2UN 

0.0901 

30* 

— 

3UN 

0.0687 

45* 

~ 

lUR 

0.0471 

45* 

0.0403 

2UR 

0.0470 

45* 

OXMOO 

3UR 

0.0376 

45* 

0.0305 

IBN 

0.1390 

30* 

2BN 

0.1390 

30* 

~ 

3BN 

0.0992 

30* 

“ 

IBR 

0.0669 

30* 

0.0618 

2BR 

0.0669 

30* 

0.0618 

3BR 

0.0514 

30* 

0.0441 

Cue  Dewriiiliaa; 


1  -  Open  hole,  2  •  bleifeieooe  Ht  Futenee,  3  -  Bonded  been, 

U .  UnUziel  Loed,  B  •  Biaxiil  Uwd.  N  -  No  Rq>eir,  R  -  Repeiicd. 


•  • 
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Tabic  19:  Resultt  of  lUiic  compression  tests:  Single  impact  case 


SptemtH  Impacior 

Absorbed 

Damaged 

Unrepaired  (U) 

Faiiure 

Failure 

Predicted 

mmbtr 

diameter 

energy 

area 

Reqtured  (R) 

load 

strain 

strain 

(mm) 

(J) 

(mn^) 

IkN) 

<m 

Ote; 

1 

19.8 

7J9 

453 

U 

-191 

-4503 

2 

19.8 

7J5 

453 

U 

-214 

-4680 

— 

3 

19.8 

7.96 

453 

u 

-213 

-4826 

— 

4 

19.8 

5.60 

479 

R 

-238 

-4993* 

-6164 

5 

19.8 

8.21 

428 

R 

-238 

-4993* 

-6164 

6 

19.8 

7.84 

448 

R 

-289 

-6699 

-6164 

7 

19.8 

7.88 

458 

R 

-289 

-6699 

•6164 

8 

30.0 

IM 

733 

U 

-168 

-4097 

9 

30.0 

6.93 

718 

U 

-192 

-4375 

— 

10 

30.0 

7.88 

665 

u 

-173 

-4305 

— 

11 

30.0 

7.54 

761 

u 

-1% 

-4350 

— 

12 

30.0 

7.74 

761 

u 

-197 

-4274 

— 

13 

30.0 

7.67 

761 

u 

-178 

-4025 

— 

14 

30.0 

7.66 

800 

R 

-233 

-5293 

-5594 

15 

30.0 

7.35 

800 

R 

-233 

-5293 

-5594 

16 

30.0 

7.45 

739 

R 

-227 

-5554 

-5594 

17 

30.0 

7.64 

704 

R 

-227 

-5554 

-5594 

18 

39.7 

6.25 

1252 

U 

-178 

-4061 

19 

39.7 

5.76 

1252 

u 

-204 

-4445 

.. 

20 

39.7 

5.87 

1252 

u 

-181 

-4090 

.. 

21 

39.7 

7.83 

1385 

R 

-238 

-5337 

-5542 

22 

39.7 

7.58 

1212 

R 

-238 

-5337 

-5542 

23 

39.7 

7.70 

1290 

R 

-222 

-5099 

-5542 

24 

39.7 

7.59 

1120 

R 

-333 

-5099 

-5542 

S|ifjdnKin  exhibited  extensive  bending  prior  to  fxihite  (anti-biickliiig  rig  wu  distorted) 


•  • 
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TaUc  11:  Strain  gauga  itaults  for  the  damagnd  and  undamaged  skins  on  the  CREDP 
hosizontal  stabilaior.  CiHnprcsaion  case  measurements  in  lu!) 


RBM(kNm) 

0 

-9 

-18 

-25 

-31 

-25 

-18 

-9 

0 

Torqut  (kNm) 

0 

-6 

-13 

-18 

-22 

-18 

-13 

-6 

Gauge  D1 

+45 

-1 

-131 

-258 

-361 

-453 

-369 

-269 

-144 

-1 

Gauge  D1 

0 

-1 

-263 

-516 

-717 

-899 

-732 

-531 

-283 

-1 

Gauge  D1 

-45 

0 

-49 

-95 

-132 

-164 

-132 

-96 

-51 

-1 

Gauge  D2 

+45 

-1 

-64 

-125 

-175 

-219 

-176 

-127 

69 

-1 

Gauge  D2 

0 

0 

-119 

-232 

-321 

-400 

-327 

-240 

-129 

-1 

Gauge  D2 

-45 

-1 

-12 

-24 

-32 

-39 

-33 

-30 

-17 

-1 

Gauge  D3 

+45 

0 

•182 

-356 

-496 

-622 

-505 

-366 

-194 

-1 

Gauge  D3 

0 

- 

423 

-827 

-1151 

-1437 

-1156 

-831 

-434 

-2 

Gauge  D3 

-45 

0 

-95 

-186 

-259 

-325 

-262 

-191 

-101 

-2 

Gauge  D4 

+45 

-1 

-126 

-246 

-344 

-431 

-351 

-256 

-138 

-1 

Gauge  D4 

0 

-1 

-261 

-508 

-705 

-878 

-717 

-522 

-281 

-1 

Gauge  D4 

-45 

-1 

-25 

-47 

-63 

-78 

-645 

-48 

-27 

-1 

Gauge  Ul 

+45 

0 

119 

232 

322 

401 

328 

238 

127 

-1 

Gauge  UI 

0 

-1 

260 

507 

705 

877 

714 

518 

275 

-1 

Gauge  Ul 

-45 

0 

60 

119 

165 

204 

164 

118 

60 

-1 

Gauge  U2 

+45 

-1 

113 

219 

304 

376 

308 

224 

119 

-1 

Gauge  U2 

0 

-1 

161 

316 

440 

550 

446 

320 

168 

-1 

Gauge  U2 

-45 

-1 

-16 

-27 

-36 

-42 

-39 

-37 

-24 

-1 

Gauge  U3 

+45 

-1 

171 

333 

463 

573 

467 

339 

180 

-1 

Gauge  U3 

0 

-1 

405 

790 

1097 

1357 

1107 

804 

429 

-1 

Gauge  U3 

-45 

-1 

112 

218 

302 

373 

304 

220 

117 

-1 

Gauge  U4 

+45 

-1 

93 

182 

252 

311 

252 

181 

95 

-1 

Gauge  U4 

0 

-1 

257 

502 

702 

872 

708 

513 

272 

0 

Gauge  U4 

-45 

0 

-3 

-6 

-8 

-10 

-8 

-6 

-3 

0 

U 

1 

370 

734 

1038 

1285 

1067 

788 

438 

0 

D 

0 

-363 

-719 

-1016 

•1260 

-1041 

-754 

^3 

0 

•  •! 
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Table  12:  Laminate  properties  for  the  horizontal  stabilaior  skin  and  repair  doubln- 


Exx 

Err 

Gxr 

vxr 

vyx 

(GPa) 

(GPa) 

(GPa) 

Skin 

66.7 

32.5 

■Ql 

0.389 

0.190 

Repair 

57.4 

32.9 

0.443 

0.253 

Table  13:  Summary  of  strain  results  from  finite  element  analysis  at  DLBM 


Strain  gauge 
number 

Unrepaired 

pe 

Repaired 

M* 

% 

Reduction 

1 

4,953 

3,557 

28 

2 

4,621 

3399 

26 

3 

5,210 

3,658 

30 

4 

4,863 

3322 

28 

LVDT 

5^84 

3,664 

31 

Figure  1. 


Plane  view  of  the  delaminaiion. 
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Z/'/y. 

Plane  of  symmetry  • 

(b) 


Figure  2.  Cross-sectional  view  of  the  delamination.  (a)  Specimen  A,  (b)  Specimen  B. 
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Figure  3.  Typical  load/deflection  curve. 


Figure  4.  Variation  of  G  with  radius. 


Figure  6.  Mesh  geomciry. 
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Figure  7.  Plot  of  analytical  and  numerical  values  of  G  with  increasing  load. 


Analytical  solution:  - 

Displacamant  sxtrapolation  method: 
(4-node  mesh)  •,  and  (8-node  mesh)  □. 
P  =  200kN 


Figure  8.  Plot  of  analytical  and  numerical  values  of  G  with  radius. 
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Figure  9.  Dculils  of  finilc  element  model. 


Figure  10.  (a)  C -.scan  of  initial  damage  from  |15|.  (b)  C-scan  of  Final 
dclamination  growth  after  20x  10^  cycles  1 1.*)!. 
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Figure  12.  Influence  of  shear  moduli  on  the  r*-iniegral. 
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Figure  14.  Distribution  of  3'*pi.«io/^ej»«ic  "'****  increasing  plasticity. 
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